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A  second  thrust  of  our  experimental  program  has  been  a  study  of  electron 
beam  induced  optical  emission  in  helium.  In  the  course  of  this  research  we 
discovered  a  new  effect  which  we  published  in  Applied  Physics  Letters. 
Relativist!^  electron  beam  interactions  with  very  small  ratios  o£  n^tr^geg  t^>. 
helium  (10  -10  )  have  been  found  to  produce  extremely  large  (B^> 

intensities  at  391.4  and  427.8  nm,  compared  to  line  intensities  originating® 
from  helium.  These  results  occurred  in  the  total  pressure  regime  of  0.1-500 
Torr.  The  pressure  scaling  results  discovered  are  inconsistent  with  previously 
proposed  kinetic  mechanisms  for  the  laser  pumped  by  helium.  With  a  simple 

model  of  the  chemical  kinetics,  we  have  shown  that  this  effect  is  due  to  the 
collisional  transfer  of  energy  between  excited  states  of  helium  atoms  and  the 
ground  state  of  N„  . 

The  main  goal  of  the  theoretical  investigation  is  to  formulate  a  model 
which  will  describe  the  motion  of  the  carbon  plasma  and  predict  the  spectrum 
of  emitted  light.  The  predicted  spectra  may  then  be  directly  compared  to 
spectra  which  are  experimentally  measured.  The  model  being  currently  developey 
is  based  on  previous  work  done  by  J.  W.  Thornhill  under  AFOSR  sponsorship  at 
Michigan.  Three  main  areas  must  be  addressed:  hydrodynamic  description  of  the 
plasma,  an  ionization  dynamics  model,  and  the  energy  deposition  of  the  inciden 
particle  beam.  Due  to  the  complexity  of  the  problem  numerical  methods  are 
used  to  obtain  results. 

•The  main  hydrocjynamic  assumption  is  that  the  plasma  can  be  described 
using  single  fluid  and  temperature  equations./  This  places  a  constraint  on  the 
time  scale  of  the  problem.  Namely  that  the  electron-ion  thermal izat ion  time 
should  always  be  smaller  than  the  hydrodynamic  time.  This  constraint  is  not 
serious  but  should  be  kept  in  mind.  The  model  is  developed  in  one  dimensions  1 y 
Lagrangian  coordinates.  This  simplifies  the  conservation  equations  and  allows 
for  simulation  of  situations  where  extensive  compression  or  expansion  occurs. 

•  The  Co  1 1  i  s  ional- Rad  ia  t  i  ve  Equilibrium  (CRE)  model  is  used  to  describe  the 
ionization  dynamics  of  the  system.  This  model  consists  of  a  set  of  rate 
equations,  one  for  each  quantum  state  included,  which  describe  the  effect  of 
atomic  processes  on  the  population  and  depopulation  of  a  given  state.  The 
model  used  is  an  equilibrium  model  since  no  time  dependence  is  included  in  the 
rate  equat ions . 

'‘The  final  component  of  our  theory  is  the  energy  depositon  model.  The 
pproach  is  simplistic  at  the  present.  The  stopping  power  of  an  incident  \']t_ 

relativistic  electron  is  integrated  over  the  plasma.  This  is  tnen  multi¬ 
plied  by  the  cur  rent  dens  ity  of  the  beam  at  each  point  in  the  plasma  to  obtain 
the  desired  result.  'JThe  stopping  power  is  based  on  a  formalism  for  r  e  la  t  iv  i  s- 
tic  electrons  interacting  with  a  plasma  which  has  both  free  and  bound  elec-  cl' 
trons.  A  collision  term  was  found  using  Fokker- Planck  formalism.  Assuming 
the  incident  beam  can  be  treated  as  a  delta  function  in  momentum  space,  the  t 
appropriate  moment  of  the  collision  term  yields  the  stopping  power. 
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We  have  performed  extensive  experiments  to  investigate  the 
spectroscopic  emission  from  electron  beam  driven  carbon  anode  plasmas. 
Diode  closure  results  in  three  phases  of  beam-plasma  evolution  which  are 
summarized  below: 

Stage  I:  Deposited  electron  dose  is  smaller  than  that  required  for 
anode  plasma  formation.  No  light  is  detected  in  this  phase. 

Stage  II:  Deposited  electron  dose  is  large  enough  for  anode  plasma 
formation  before  the  diode  shorts.  Spectroscopic  data  for 
noncrowbarred  shots  (dose  »  500-800  J/g)  indicate  the  presence  of 
CII  and  CIII  with  a  high  level  of  background  radiation. 

Stage  III:  Cathode  plasma  begins  shorting  A-K  gap  or  anode  and 
cathode  plasmas  meet.  For  crowbarred  shots  (dose  >  300  J/g) 
molecular  emission  is  observed  from  C2  and  CH.  Atomic  emission  from 
CII  and  H.  Noncrowbarred  shots  exhibit  intense  atomic  emission;  at 
early  times  CIII  and  CII  are  the  primary  components,  however 
emission  from  CIV  has  also  been  observed.  Electron  temperatures 
from  Boltzmann  plots  are  2-3  eV.  At  later  times  the  emission  is 
primarily  from  Cl  and  CII  with  Te  *  2-3  eV.  The  electron  density  is 
found  to  be  Ne=1018  cm'3  from  the  Saha  equation. 

A  second  thrust  of  our  experimental  program  has  been  a  study  of 
electron  beam  induced  optical  emission  in  helium.  In  the  course  of  this 
research  we  discovered  a  new  effect  which  we  published  in  Applied 
Physics  Letters.  Relativistic  electron  beam  interactions  with  very  small 
ratios  of  nitrogen  to  helium  (10'1-10'4)  have  been  found  to  produce 
extremely  large  N2+(B2Iu+-X2Ig+)  intensities  at  391.4  and  427.8  nm, 
compared  to  line  intensities  originating  from  helium.  These  results 
occurred  in  the  total  pressure  regime  of  0.1-500  Torr.  The  pressure 
scaling  results  discovered  are  inconsistent  with  previously  proposed 
kinetic  mechanisms  for  the  N2+  laser  pumped  by  helium.  With  a  simple 
model  of  the  chemical  kinetics,  we  have  shown  that  this  effect  is  due  to 


the  collisional  transfer  of  energy  between  excited  states  of  helium  atoms 
and  the  ground  state  of  N2+- 

The  main  goal  of  the  theoretical  investigation  is  to  formulate  a  model 
which  will  describe  the  motion  of  the  carbon  plasma  and  predict  the 
spectrum  of  emitted  light.  The  predicted  spectra  may  then  be  directly 
compared  to  spectra  which  are  experimentally  measured.  The  model  being 
currently  developed  is  based  on  previous  work  done  by  J.W.  Thornhill  under 
AFOSR  sponsorship  at  Michigan.  Three  main  areas  must  be  addressed: 
hydrodynamic  description  of  the  plasma,  an  ionization  dynamics  model,  and 
the  energy  deposition  of  the  incident  particle  beam.  Due  to  the  complexity 
of  the  problem  numerical  methods  are  used  to  obtain  results. 

The  main  hydrodynamic  assumption  is  that  the  plasma  can  be 
described  using  single  fluid  and  temperature  equations.  This  places  a 
constraint  on  the  time  scale  of  the  problem.  Namely  that  the  electron-ion 
thermalization  time  should  always  be  smaller  than  the  hydrodynamic  time. 
This  constraint  is  not  serious  but  should  be  kept  in  mind.  The  model  is 
developed  in  one  dimensional  Lagrangian  coordinates.  This  simplifies  the 
conservation  equations  and  allows  for  simulation  of  situations  where 
extensive  compression  or  expansion  occurs. 

The  Collisonal-Radiative  Equilibrium  (CRE)  model  is  used  to 
describe  the  ionization  dynamics  of  the  system.  This  model  consists  of  a 
set  of  rate  equations,  one  for  each  quantum  state  included,  which  describe 
the  effect  of  atomic  processes  on  the  population  and  depopulation  of  a 
given  state.  The  model  used  is  an  equilibrium  model  since  no  time 
dependence  is  included  in  the  rate  equations. 

The  final  component  of  our  theory  is  the  energy  deposition  model. 
The  approach  is  simplistic  at  the  present.  The  stopping  power  of  a 
incident  relativistic  electron  is  integrated  over  the  plasma.  This  is  then 
multiplied  by  the  current  density  of  the  beam  at  each  point  in  the  plasma 
to  obtain  the  desired  result.  The  stopping  power  is  based  on  a  formalism 
for  relativistic  electrons  interacting  with  a  plasma  which  has  both  free 
and  bound  electrons.  A  collision  term  was  found  using  Fokker-Planck 
formalism.  Assuming  the  incident  beam  can  be  treated  as  a  delta  function 
in  momentum  space,  the  appropriate  moment  of  the  collision  term  yields 
the  stopping  power. 


1.0)  Spectroscopy  of  Electron  Beam  Induced  and  Driven  Carbon  Plasmas 


1.1)  Experimental  Configuration 

These  experiments  are  being  performed  on  MELBA  (Michigan  Electron 
Long  Beam  Accelerator)  at  the  University  of  Michigan.  MELBA  is  driven  by  a 
state  of  the  art,  long  pulse  Marx  bank  with  a  special  Abramyan, 
reverse-charged  capacitor  stage  for  voltage  compensation.  Typical  MELBA 
parameters  that  have  been  obtained  are  peak  voltages  <  1  MV,  currents  <  30 
kA,  and  pulselengths  =  0.3  -  4.0  ps.  In  this  section  the  experimental 
equipment  and  configuration  are  described.  In  addition  some  improvements 
are  noted  for  future  experiments  where  appropriate. 

The  principal  idea  in  the  Abramyan  concept  is  illustrated  in  Figure  1.1. 
The  voltage  Vk(t)  is  the  reverse-charged  stage  voltage.  This  is  positive 
initially,  with  a  maximum  of  about  60  kV,  subtracting  from  the  main  Marx 
voltage  Vg(t).  Note  Vg(t)  droops  due  to  the  natural  RC  decay  of  the  Marx. 
This  decay  is  exacerbated  with  a  cold  cathode  diode  load  due  to  collapsing 
diode  impedance.  Later  in  the  pulse  the  reverse  charged  stage  rings 
negative,  adding  on  to  the  main  Marx  voltage.  The  combination  of  initial 
subtraction  and  later  addition  of  the  Abramyan  voltage  is  intended  to  give 
the  flat  total  voltage  pulse  Vg(t)  +  Vk(t).  The  design  goals  were  to 
compensate  for  an  impedance  collapse  of  a  diode  from  300  Q  to  50  D  in 
1.075  ps,  and  giving  a  voltage  flat  to  within  ±10%.  Figure  1.2  shows 
comparison  of  experiment  with  theory  for  compensation  of  the  RC  decay 
into  a  127  Q  fixed  resistive  load.  Compensation  has  also  been 
demonstrated  with  large  gap  diodes.-*-1 

The  MELBA  circuit  with  diode  load  is  shown  in  Figure  1.3.  It  displays 
several  important  machine  features  in  addition  to  the  Abramyan  stage, 
here  represented  by  the  elements  RA,  CA,  VA,  LA.  The  crowbar  circuit,  in 
parallel  with  the  diode  load,  is  designed  to  switch  the  voltage  away  from 
the  diode  by  shunting  the  majority  of  the  current  to  ground.  The  crowbar 
gap,  a  high  pressure  SF6  switch,  can  be  command  triggered  with  different 
time  delays  in  order  to  obtain  different  pulselengths.  The  diode  load  is 
represented  by  the  stalk  inductance  Ls  =  400  nh,  and  a  time  changing  diode 


MELBA  OUTPUT  VOLTAGE  INTO  1 27  OHM  RESISTIVE  LOAD 


Figure  1.2 
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impedance  ZD(t)  with  both  inductive  and  resistive  components.  The  other 
important  feature  is  the  presence  of  a  large  series  resistance  Rseries  = 

Q,  intended  to  take  the  full  generator  energy  in  a  crowbar  or  in  the  event 
of  a  diode  short  to  prevent  Marx  damage.  This  resistor  plays  an  important 
role  in  determining  load  dynamics  when  RD  <  RSeries-  a  case  which  occurs 
for  smaller  gaps  and  late  in  time  after  gap  closure  has  occurred. 

The  typical  stored  energy  in  the  Marx  bank  is  25  kJ.  In  experiments 
with  large  diode  gaps  (8-10  cm)  energies  of  10  kJ  have  been  delivered  to 
the  anode  at  low  current  densities  (<20  A/cm2).  The  experiments 
discussed  in  this  report  were  performed  at  intermediate  gaps  (3-4  cm) 
and  have  delivered  up  to  7  kJ  at  current  densities  one  to  two  orders  of 
magnitude  larger  (>100  A/cm2).  The  present  study  indicates  that  the 
critical  parameter  is  not  the  total  delivered  energy  but  the  energy  density 
or  dose,  D  (J/gr),  which  requires  knowledge  of  the  current  density.  In  the 
large  gap  experiments  D  =  20  J/gr.  The  dose  for  the  intermediate  gap 
experiments  was  increased  by  one  to  two  orders  of  magnitude,  resulting  in 
more  intense  atomic  optical  emission. 

The  most  recent  improvements  in  the  experimental  configuration  were: 

i)  higher  electron  current  densities  by  using  shorter  gaps, 

ii)  exclusive  use  of  noncrowbarred  shots, 

iii)  extended  anode  to  view  plasma  immediately  on  the  face, 

iv)  and  different  cathode  and  anode  materials  to  distinguish 
cathode  plasma  from  anode  plasma. 

Further  improvements  in  current  density  from  gap  reduction  can  still 
be  made.  In  addition,  we  have  designed  pulsed  magnetic  field  coils  for  the 
diode  which  have  been  built  by  an  outside  contractor.  We  are  installing 
them  on  the  diode  chamber  .  The  maximum  field  on  axis  expected  will  be 
10  KG.  The  externa!  field  will  result  in  increases  in  current  density. 

The  principle  reason  that  noncrowbarred  shots  are  so  desirable  is  that 
the  current  density  increases  very  quickly  late  in  time  as  the  plasma  nears 
closure.  This  greatly  enhances  the  anode  dose.  The  difference  between 


crowbarred  and  noncrowbarred  shots  is  an  increase  in  total  dose  by  a 
factor  of  two  to  four  by  the  end  of  the  pulse.  The  use  of  different  anode 
and  cathode  materials  in  order  to  distinguish  the  plasma  source  is  a 
significant  advantage.  The  cathode  plasma  components  are  pump  oil  and 
hydrocarbon  impurities  C2,  CH,  H,  as  well  as  neutral  and  singly  ionized 
carbon  and  cathode  metal  components.  The  carbon  components  from  the 
cathode  plasma  do  not  interfere  with  distinguishing  carbon  produced  from 
the  anode  as  they  are  much  cooler  and  contain  significant  molecular 
constituents. 

The  diode  configuration  is  shown  in  Figure  1.4.  The  cathodes  used  were 
made  of  either  brass  or  aluminum,  with  a  small  cross  sectional  area  and 
had  a  blunt  nose.  One  method  of  making  further  improvements  in  the 
incident  current  density  is  to  make  the  cathode  have  a  sharp  point.  Figure 
1.4  also  shows  the  extended  anode. 

The  method  of  light  detection  and  measurement  is  shown  in  Figure  1.5. 
The  initial  collecting  lens  is  intended  to  collimate  the  source  light  which 
is  transmitted  to  the  optical  multichannel  analyzer  (OMA)  via  mirrors  and 
a  final  focusing  lens.  Without  taking  into  account  the  finite  optical  depth 
of  the  lenses,  the  system  has  a  spatial  resolution  of  about  100  pm,  the 
size  of  the  final  slit  (MT  =  1).  In  experiments  thus  far  the  light  has  never 
been  perfectly  collimated.  The  acrylic  windows  we  are  using  in  the  current 
experiments  have  a  cut  off  for  transmission  below  about  3800  A.  Future 
experiments  at  shorter  wavelengths  will  require  quartz  optics  and 
windows  and  an  evacuated  optical  path.  Some  of  the  quartz  lenses  have 
been  already  purchased.  The  acrylic  windows  become  clouded  with  a  film 
by  the  carbon  vapor  and  cannot  be  cleaned  with  conventional  methods.  The 
quartz  windows  can  be  recovered  to  new  optical  quality  by  baking  them  at 
450  °C  for  24  hours. 

The  OMA  equipment  is  housed  in  a  lead  faced  screen  room  to  prevent 
damage  to  the  sensitive  electronics  by  hard  x-rays.  The  OMA  consists  of  a 
0.275  meter  Jarell  Ash  spectrograph  and  a  Tracor  Northern  1024  channel 
gated  intensified  diode  array.  Plans  are  being  made  to  replace  the  0.275 
meter  spectrograph  with  a  1.0  meter  version  (supported  by  AFOSR),  greatly 
increasing  our  resolution  capabilities.  For  example,  typical  best  resolution 
for  the  0.275  meter  spectrograph  is  about  10  A.  This  is  expected  to  be 
reduced  to  0.2  A  with  the  1.0  meter  model. 
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1 .2)  Experimental  Modeling 


1 ,2.a)  Diode  Analysis 

One  of  the  principle  features  of  this  beam-plasma  experiment  different 
from  others  is  that  the  interaction  occurs  in  the  diode  region.  Plasma 
formed  in  the  diode  affects  the  electron  beam  coupling  as  a  result  of  the 
diode  dynamics  and  the  diode-generator  coupling.  For  example,  the 
formation  of  an  anode  plasma  in  an  electron  diode  can  result  in  an 
increased  electron  current  due  to  space  charge  neutralization.  It  may  also 
short  out  the  diode  faster  if  it  has  a  finite  velocity.  Another  possibility  is 
that  it  may  increase  the  demands  on  the  generator  or  the  cathode  plasma 
more  quickly  or  past  the  point  where  they  can  be  supplied,  resulting  in  an 
opening  of  the  gap,  rather  than  continued  impedance  collapse.  Thus  diode 
physics  is  an  important  part  of  this  study  in  determining  details  about  the 
beam  plasma  coupling. 

An  idealized,  pictorial  example  of  the  diode  time  development  is  shown 
in  Figure  1.6.  The  position  of  the  OMA  gate  relative  to  these  possible  light 
producing  events  in  the  diode  determines  the  source  of  the  detected  light. 
Thus  diode  physics  is  also  important  in  order  to  determine  what  is 
happening  in  the  diode  when  the  OMA  is  gated.  An  extended  analysis  and 
modeling  o‘  the  voltage  and  current  is  required.  For  the  purposes  of  this 
discussion  and  modeling,  the  evolution  of  the  diode  proceeds  through  three 
stages,  defined  in  Figure  1.7  and  as  follows. 

STAGE  I  -  Initial  stages  of  cathode  plasma  formation  and  expansion.  The 
anode  dose  is  below  that  required  for  anode  plasma  formation. 
No  light  would  be  detected. 

STAGE  II  -  In  this  stage  the  dose  exceeds  some  minimum  required  dose 
for  anode  plasma  formation  before  the  cathode  plasma  has 
completely  crossed  the  gap  and  shorted  the  diode.  In  response 
to  the  presence  of  anode  ions,  we  expect  a  bipolar  current.  The 
anode  plasma  may  or  may  not  close.  Some  shots  may  not  have 
this  stage.  In  electron  diodes,  carbon  anode  plasmas  are 


IDEALIZED  DIODE  TIME  DEVELOPMENT 


THE  POSITION  OF  THE  OMA  GATE  RELATIVE  TO  THESE  POSSIBLE  LIGHT 
PRODUCING  EVENTS  DETERMINES  THE  SOURCE  OF  THE  DETECTED  LIGHT. 
AN  EXTENDED  ANALYSIS  IS  REQUIRED  IN  ORDER  TO  ASSESS  WHICH  ARE 
OCCURING  DURING  THE  GATE. 

STAGE  I  -  DOSE  SMALLER  THAN  THAT  REQUIRED  FOR  ANODE  PLASMA  FORMATION 


STAGE  II  -  DOSE  LARGE  ENOUGH  TO  FORM  ANODE  PLASMA  BEFORE  SHORTING 
(SOME  SHOTS  MAY  NOT  HAVE  THIS  STAGE) 


STAGE  III 

CATHODE  PLASMA  BEGINS  SHORTING  ANODE  AND  CATHODE  PLASMAS 

MEET  AND  SHORT  GAP 


expected  to  be  formed  between  350  and  600  J/gr.  These 
plasmas  are  formed  from  ionized  carbon  vapor  from  a  heated 
outgassing  carbon  anode.  The  vaporization  dose  for  direct 
carbon  sublimation  is  =  56  kJ/gr  at  a  temperature  of 
3600  °C. 

STAGE  III-  If  no  anode  plasma  has  formed,  then  the  cathode  plasma  starts 
to  short  out  the  gap.  As  the  plasma  gradually  fills  the  whole 
diode  and  touches  the  entire  anode  surface,  the  emission 
surface  for  space  charge  emitted  electrons  is  gradually  shut 
off  and  the  diode  acts  in  more  of  a  glow  discharge  or  vacuum 
arc  mode.  This  is  a  transition  region  where  both  an  electron 
beam  and  an  arc  exist.  For  shots  which  have  had  a  stage  II,  the 
anode  and  cathode  plasmas  meet  and  the  same  type  of  gradual 
transition  takes  place.  We  distinguish  between  early,  middle 
and  late  divisions  of  stage  III.  Early  stage  III  contains  the 
transition  region,  thus  is  a  mixture  of  beam  and  arc. 

What  follows  is  a  discussion  of  the  diode  physics  necessary  to 
determine  the  events  in  the  diode  and  the  stage  for  a  particular  shot.  It 
includes  a  discussion  of  diode  closure  analysis,  a  simple  diode  current 
model  and  an  anode  dose  model.  The  diode  closure  analysis  yields  the 
cathode  closure  velocity  u,  (an  input  to  the  current  model  and  the  dose 
model),  and  effective  diode  area,  (an  input  to  the  current  model).  The 
closure  velocity  analysis  gives  the  expected  cathode  plasma  diode  shorting 
time  Tshort.  The  position  and  width  of  the  OMA  gate  relative  to  this  time 
is  important  for  two  reasons.  First  it  gives  an  estimate  of  the  time  to 
expect  the  transition  to  the  arc  mode,  and  secondly  when  to  expect  cathode 
plasma  components  to  become  visible.  We  must  keep  in  mind  that  this 
velocity  is  only  an  apparent,  average  uniform  velocity  for  the  cathode 
plasma.  The  actual  plasma  probably  has  faster  and  slower  components  and 
is  nonuniform.  Use  of  this  velocity  is  valid  for  comparison  since  the  same 
method  is  applied  for  each  shot. 

The  ability  to  model  the  diode  current  is  important  in  order  to  verify 
that  the  diode  is  behaving  in  a  mode  we  understand,  and  to  ascertain  when 
bipolar  effects  may  come  into  play.  This  would  give  some  confirming 


evidence  for  the  presence  of  an  anode  plasma.  Finally,  the  dose  model,  here 
only  a  crude  estimate  of  the  incident  electron  energy  density,  again  gives 
some  indication  of  whether  to  expect  an  anode  plasma  during  the  OMA  gate, 
and  to  tentatively  decide  what  stage  the  shot  is  in.  The  dose  model  is  more 
certain  in  the  portion  of  the  pulse  that  can  be  modeled  well  by  space 
charge.  Again  the  uncertainty  in  the  dose  model  is  unimportant  as  all  shots 
are  treated  identically.  An  accurate  determination  of  the  dose  would 
require  an  actual  measurement  of  current  density,  rather  than  current. 
Following  the  theory  discussion  is  a  sample  application  to  modeling 
idealized  data. 

i)  Diode  Closure 

The  Child  -Langmuir  law  for  a  space  charge  limited  infinite  planar 
diode  is  given  by: 

ICL  =  2.336  XI 0‘6  n  (KRr)2V3/2  (1.1) 

"(d0-uc(T-TSTART))^ 

In  this  equation,  d0  is  the  initial  gap  spacing,  uc  is  the  cathode  plasma 
closure  velocity,  Rc  is  the  cathode  radius,  Tstart  is  the  start  of  the  Child 
Langmuir  behavior  region  or  the  diode  opening  time,  and  K  is  a  constant 
(*1)  to  account  for  relativistic,  nonplanar,  and  ion  effects. 

To  perform  the  closure  analysis  solve  Equation  1.1  for  the  denominator: 

deff(t)  =  d0  -  uc(T  -  Tstart)  -  2.709  X  10-3Reff(1/P)1/2  (1.2) 

In  this  equation  P  «  lexp/Vexp3/2,  and  Ref{  =  KRC.  Now  perform  a  linear  best 
fit  to  the  quantity  (l/P)"^  determined  from  experimental  data.  The  fit 
should  be  performed  between  the  times  Tstart  and  TEND.  These  times  are 
chosen  so  that  the  fit  has  a  high  correlation  coefficient  which  signifies 
that  the  space  charge  flow  model  is  expected  to  give  a  good  match  to  the 
data.  The  best  fit  yields  the  equation  (1/P)fjt1/2  =  A  -  B(T  -  Tstart)  where 
A  and  B  are  the  intercept  and  slope  respectively  of  the  best  fit  line  to 
points  between  Tstart  and  TEND.  The  (l/P)1^  best  fit  analysis  is  depicted 


in  Figure  1.7.  Substitute  this  fit  for  (1/P)1/2  in  Equation  1.2  to  obtain  the 
best  fit  representation  for  deff(t).  Now  define  TSHOrt  -  d0/uc  +  Tstart.  At 
T  =  Tshort  deff  is  defined  to  be  zero.  Use  this  property  in  the  best  fit 
representation  for  deff(t)  above  and  solve  for  uc: 

Uc  =  d0B/A  (1.3) 

Other  estimates  of  the  shorting  time,  for  example  deff  is  zero  at  TEND  will 
yield  the  same  estimate  for  uc.  The  quantity  Re{{  can  also  be  determined 
from  this  same  best  fit  analysis.  Assuming  that  at  TStart-  deff  is 
physically  equal  to  the  initial  gap  spacing  d0  yields: 

Refffit  =  KRc  =  do  0-4) 

2.709x1 0'3A 

ii)  Current  Model 

The  diode  closure  analysis  yields  both  uc  and  Refffit-  Thus  a  model  for 
the  diode  perveance  during  the  space  charge  portion  of  the  pulse  is 
completely  determined  by  the  linear  best  fit  to  (1/P)1^. 

*exp(*)  =  Pexp(t)  =  Pmodel(t)  =  2.336x10 6  n  Reff2  (1.5) 

VexpW372  (do  ‘  uc(T '  tstart))2 

The  currents  predicted  by  this  model  match  the  experimental  currents 
well  until  near  TEND.  Recall  that  TEND  was  chosen  to  give  the  best  fit  with 
the  highest  correlation  coeficient.  The  match  of  the  current  model  to 
measured  current,  the  real  indicator  of  pure  space  charge  behavior,  is  now 
used  to  fine  tune  the  position  of  TEND.  TEND  is  adjusted  so  that  for  T  > 
Tend  the  unipolar  space  charge  model  does  not  match  the  experimental 
data.  Again  this  procedure  is  illustrated  in  Figure  1 .7. 

Another  informative  rearrangement  of  the  space  charge  equation 
involves  solving  for  the  Reff  required  in  order  to  match  the  instantaneous 
experimental  current: 


Reff(*)  *  369.1 4[d0  -  uc(T  -  Tgy^y)]  [Pgxp^)] 


(1.6) 


1/2 


The  best  fit  Reff  oscillates  around  the  instantaneous  Reff(t)  during  the 
portion  of  the  pulse  where  the  space  charge  model  is  a  good  match  (as 

expected).  After  TEND  .  Reff(t)  decreases  towards  zero,  a  decrease  in  the 

effective  diode  perveance,  corresponding  to  a  change  in  diode  operating 
mode,  possibly  representing  the  transition  to  an  arc.  The  diode  or  load 
resistance  continues  to  drop  during  the  transition.  At  this  point  the 

resistance  of  the  load  is  only  a  few  ohms.  In  the  limit  of  zero  load 

resistance,  the  series  resistor  limits  the  generator  current  to  s  V/4  <  25 
KA  for  voltages  late  in  the  pulse.  This  current  is  approximately  the  limit 
on  late  time  current  measured.  The  shot  to  shot  variation  of  Reff  and  u  is 
large.  Physically  this  represents  a  changing  effective  diode  load.  Thus  the 
reproducability  of  the  shots  is  poor.  This  is  another  reason  why 
comparisons  require  an  extended  analysis. 

iii)  Dose  Analysis 

The  dose  is  calculated  by  performing  the  integral: 


D(t) 


Tend 


r 

J  o 


J(t)V(t)  dt  + 
R[V(t)] 


Tf 

f  mm  dt 


J  TendAR[V(1)] 


(1.7) 


where  J(t)  is  given  by  Child-Langmuir  scaling  with  best  fit  closure 
velocity  (Equation  1.1  divided  by  n(KRc)2),  R[V(t)]  is  the  CDSA  range  in 
gm/cm2  fit  as  a  function  of  voltage,  Tp  is  the  limit  of  recorded  data,  and 
A  is  the  plasma  cross  sectional  area  =  nRc2.  This  formulation  neglects 
conduction  effects  in  the  affected  electron  maximum  range  thick  anode 
layer,  a  valid  assumption  for  the  short  timescale  of  the  pulselength.  The 
first  integral  contains  the  contribution  from  space  charge  flow  in  the 
diode.  The  second  is  a  conservative  estimate  of  the  dose  from  the  period  of 
the  pulse  were  space  charge  flow  for  the  current  is  not  a  good  model 


(after  Tend)-  This  Period  contains  the  transition  region  and  the  full  arc 
mode.  The  first  part  of  the  integral  is  by  far  the  largest  contribution.  Note 
that  this  integral  requires  a  cutoff  when  Child  Langmuir  scaling  is  used 
for  J(t).  If  the  integration  is  performed  until  T  =  TShort-  the  integral  will 
diverge.  Setting  the  cutoff  at  TEND  is  a  conservative  estimate. 

An  approximate  analytical  form  of  this  analysis  predicts  that  in  the 
space  charge  phase  higher  doses  are  delivered  for  slower  closure 
velocities  for  a  given  initial  gap  spacing: 


(1.8) 


In  Equation  1.8,  D  is  in  KJ/gr,  u  in  cm/ps,  Vo,  the  peak  voltage  in  MV,  do 
the  initial  gap  spacing  in  cm,  and  C  =  Tend/Tshort  some  fraction  of  TShort 
in  order  to  make  the  result  finite.  Higher  doses  are  also  delivered  for 
smaller  initial  gaps.  Although  J(t)  increases  faster  for  higher  u,  the  limits 
of  integration  are  cutoff  sooner,  leading  to  lower  doses.  The  analysis  also 
shows  that  the  increase  in  dose  due  to  range  shortening  at  lower  voltages 
is  offset  by  the  decrease  due  to  smaller  electron  current  densities,  and 
the  dose  actually  increases  with  Vo. 

iv)  Idealized  Data 

An  idealized  version  of  typical  noncrowbarred  experimental  data  and 
analysis  is  shown  in  Figure  1.7.  This  is  an  aid  to  developing  and 
understanding  the  voltage  and  current  analysis.  There  are  several 
important  features  displayed  in  this  figure.  The  position  and  width  of  the 
OMA  gate  is  of  course  variable.  The  voltage  always  displays  the  peaking  up 
or  opening  behavior  late  in  time  as  if  the  load  impedance  were  increasing. 
In  previous  experimental  analysis  we  thought  that  the  dip  in  voltage  was 
associated  with  arcs  or  anomalous  crowbar  behavior.  The  current  modeling 
implies  that  the  peaking  behavior  is  anomalous.  The  peaking  behavior  is 
inconsistent  with  falling  diode  resistance  which  is  expected  due  to 
closure  or  because  of  anode  initiated  bipolar  effects.  Diode  inductive 
effects  in  response  to  pinching  are  not  large  enough  to  cause  such  a  large 


Dspacecharge  =  4.25  Vo  |nri/n-C11 
dou 
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voltage  increase.  Also  note  that  in  the  modeling  of  the  current,  bipolar 
effects  are  not  invoked.  However,  that  model  does  not  agree  with 
experiment,  starting  at  the  beginning  of  this  opening  behavior.  Often,  the 
dose  at  that  point  is  well  into  the  range  for  expected  anode  plasma 
formation,  which  is  consistent  with  experiments  where  light  is  measured. 
The  final  feature  to  note  is  that  the  predicted  shorting  time  is  usually 
well  after  the  end  of  the  voltage  "opening".  Some  theories  consistent  with 
these  observations  involve  the  load-generator  coupling. 

After  the  anode  plasma  is  formed  due  to  some  required  dose  being 
surpassed,  the  impedance  of  the  diode  will  drop  by  almost  a  factor  of  two, 
immediately  demanding  more  current  from  the  generator.  If  the  demand  is 
greater  than  the  generator  can  supply,  or  the  cathode  plasma  density  is  not 
large  enough  to  meet  extraction  demands,  the  cathode  and/or  anode  plasma 
emissive  surface  may  be  pushed  back  and  the  gap  opens.  This  type  of 
voltage  peaking  can  also  be  created  by  nonuniform  filling  of  the  space 
charge  potential  in  the  gap  or  by  load  mismatch  effects.  The  understanding 
of  this  behavior  is  very  important  since  it  seems  to  be  related  to  anode 
plasma  production  and  is  required  for  a  clear  interpretation  of  the  spectra. 

We  have  developed  an  interactive,  graphical  computer  code  EDIODE3, 
which  greatly  enhances  our  diode  analysis  and  understanding.  This  code  is 
based  on  the  physics  discussed  in  the  current  section.  Figure  1.8  is  a  flow 
chart  showing  the  analysis  pathway.  Figure  137  is  an  example  of  the  code 
output. 
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IDEALIZED  DATA 
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1 ,2.b)  Spectral  Analysis 


i)  Theory 

Some  of  the  principle  goals  in  plasma  emission  spectroscopy  are  to 
determine  the  plasma  constituents,  densities,  and  temperatures.  These 
properties  can  often  yield  information  about  the  state  of  the  plasma, 
whether  it  be  equilibrium  or  non-equilibrium.  This  is  an  important 
consideration  as  a  population  inversion  is  a  non-equilibrium  state.  What 
follows  is  a  brief  discussion  of  various  ways  to  model  and  assess  the 
plasma  state  and  a  description  of  a  preliminary  analysis  of  a  measured 
spectrum. 

The  state  of  a  plasma  in  complete  thermodynamic  equilibrium  (CTE) 
can  be  described  analytically.  The  principle  assumptions  for  a  rigorous 
definition  of  a  CTE  plasma  are: 

i)  Electron  and  heavy  particle  velocity  distributions  are 
Maxwellian. 


fejtee.j)  “  (me  /2rikTe  ;)3/2exp-(mfi  ;vfi  j2/2kTG  ,) 


le,r‘-u,'le,i 
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ii)  Equipartition  of  translational  energy  exists  between 
electrons  and  heavy  components. 


Te  =  Tj  =  T 


(1.10) 


iii)  The  occupation  of  bound  electronic  states  are  described  by 
the  Boltzmann  equation  (excitation  equilibrium). 


Nj  =  N  gj  exp(-Ej/kT)/P(T) 
P(T)  =  X.gj  exp(-Ej/kT) 


(1.11) 

(1.12) 


Nj/Nj  =  (gj/gj)  expREj  -  Ej)/kT] 


(1.13) 


where  Nj  and  Ni  are  the  population  densities  for  the  ion  of 
interest  of  the  upper  and  lower  state  respectively,  N  is  the 
total  density,  Ej  and  Ei  the  energies  of  the  upper  and  lower 
state  and  P  is  the  state  sum  or  partition  function.  Equation 
1.14  is  only  valid  for  an  optically  thin  medium  and  is  not 
strictly  applicable  for  CTE  plasmas.  In  that  equation,  Iij  is 
the  optical  intensity,  Xij  the  transition  wavelength,  Aji  the 
transition  probability  and  C  is  a  constant  which  depends  on 
geometry. 

iv)  The  concentration  of  densities  of  electrons,  neutral,  and 
ionized  species  are  described  by  the  Saha  equation  (ionization 
equilibrium). 

NeNz*i/Nz  *  2(2nmekT/h2)3/2[Pz+1(T)/P2(T)]exp-Ixz/kT]  (1 .15) 

NeNzti  ,k/Nz.f 6X1 021(kT)3/2[g2+1  ik/gzj]exp-[(x2+Ek-Ej)/kT]  (1.16) 

where  Nz+1  and  Nz  are  the  number  densities  for  charge  states 
z  +1  and  z  respectively  regardless  of  energy  level,  and  %z  's 
the  ionization  potential  for  charge  state  z.  In  the  second 
version,  the  Boltzmann  equation  (Eq.  1.11)  has  been  used  to 
specify  particular  excited  state  densities  Nz+i,k  and  Nz.j.  Ek 
is  the  upper  state  energy  for  the  kth  level  in  the  z+1  ion, 
while  Ej  is  the  upper  state  energy  for  the  jth  level  of  the  z 
ion. 

v)  The  radiation  field  is  described  by  a  Planck  distribution. 

Because  of  energy  loss  due  to  nonabsorption  of  radiation,  real  plasmas 
are  almost  never  in  CTE.  In  this  case  a  condition  known  as  local 
thermodynamic  equilibrium  (LTE)  may  exist  and  conditions  (i)  -  (iv)  can  be 
satisfied  pointwise  rather  then  with  a  global,  uniform  temperature. 


Condition  (v)  no  longer  applies.  In  general,  for  LTE  to  ex$ist,  collisionai 
processes  must  be  more  important  than  radiative  ones,  so  that  the 
shortfall  of  radiative  energy  does  not  matter.  Ionization  and  excitation 
must  take  place  primarily  through  collisions.  This  is  equivalent  to 
requiring  a  high  enough  electron  density.  A  more  quantitative  requirement 
for  LTE  is  that1  ^ 

Ne  >1.1  x  1014  Te1/2  (Ej  -  E{)3  (cm*3)  (1.17) 

where  AE  =  Ej  -  Ei  is  the  transition  energy  from  state  j  to  state  i  of  the 
same  ion  species  and  Te  and  AE  are  in  ev.  This  equation  can  be  rearranged  to 
find  the  principle  quantum  level  i  for  which  the  above  condition  is  satisfied. 

Ej  >  Ej  -  2.08  x  10'5Ne1/3/Te1/6  (ev)  (1.18) 

Levels  at  and  above  i  are  in  LTE,  while  those  below  are  in  a  nonequilibrium 
state  for  which  other  models  must  be  used.  This  is  known  as  a  condition  of 
partial  local  thermodynamic  equilibrium  (PLTE).  Figure  1.91-3  presents  the 
density  limits  above  which  the  ground  states  of  the  various  carbon  ions  will 
be  in  LTE.  The  conditions  1.17  and  1.18  are  the  hardest  to  satisfy  for 
transitions  involving  the  ground  state. 

These  four  remaining  conditions  can  be  used  to  define  LTE  with  the 
useful  notion  of  a  multitemperature  description.1-4’1^  The  temperatures 
associated  with  requirement  (i)  and  (ii)  are  the  species  kinetic 
temperatures.  In  a  3-species  carbon  plasma,  for  example,  one  notation  for 
the  kinetic  temperatures  would  be  Te  ,T^  J^jj.  The  conventional,  relevant 
plasma  temperature  is  of  course  the  electron  kinetic  temperature  Te.  The 
temperatures  associated  with  requirement  (iii)  are  defined  as  the 
excitation  temperatures,  in  this  example  Texcj,  TexCjj.  They  are  obtained 
from  the  Boltzmann  equation  and  the  bound  electronic  state  distribution. 
The  relevant  temperature  for  requirement  (iv)  is  defined  as  the  ionization 
temperature  TJe  and  is  obtained  from  the  Saha  equation  and  the  current 
plasma  species  concentrations.  For  the  above  example  of  a  three  component 
plasma  (electrons,  neutral  carbon-CI  and  singly  ionized  carbon-CII),  there 
are  thus  6  relevant  temperatures.  For  a  plasma  in  kinetic  equilibrium  the 
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requirements  (i)  and  (ii)  of  a  Maxwellian  particle  distribution  and 
equipartition  of  translational  energy  imply  T£  =  Td=Te  .  a  plasma  is  defined 
to  be  in  ionization  equilibrium  if  T^Tg.  Excitation  equilibrium  is  defined 
by  the  requirement  that  Texd=Texcn  =  Te.  Thus,  in  this  example,  LTE  exists 
when  Te=Tc  =  T£jj=  T*e  =  Texci  =  TexQjj=T,  i.e.  there  is  only  one  relevant 
local  thermodynamic  temperature.  Each  of  these  particular  types  of 
equilibrium  is  promoted  by  increasing  electron  density.  Ionization 
equilibrium  however  is  destroyed  by  recombination.  Both  ionization  and 
excitation  equilibrium  also  depend  on  the  various  characteristic  rates  of  the 
plasma  which  are  a  function  of  pumping  power,  plasma  hydrodynamic  time 
scales  determined  by  expansion  and  material  supply  rate,  recombination 
time  scales,  and  atomic  relaxation  time  scales.  Further  work  is  required  to 
identify  the  specific  conditions  required  for  each  type  of  equilibrium,  and  to 
evaluate  these  characteristic  rates  for  our  plasma. 

In  the  absence  of  LTE  it  is  still  possible  to  calculate  atomic  level 
populations,  but  only  if  the  relevant  collision  cross  sections  and  transition 
probabilities  are  known.  This  requires  a  great  deal  of  atomic  data  which 
may  only  be  be  known  accurately  for  H-like  and  He-like  systems.  The  next 
most  common  approximation  in  this  category  is  the  coronal  equilibrium  (CE) 
model,  which  is  appropriate  for  high  temperature,  low  density  plasmas.  For 
some  plasmas,  a  combination  of  LTE  and  CE  may  be  appropriate.  Figure 
1.1  O'*  -3  presents  the  density  limits  below  which  the  lower  states  of  the 
various  carbon  ions  are  in  CE. 

The  collisional-radiative  equilibrium  (CRE)  model  is  a  steady  state 
solution  to  the  atomic  level  rate  equations,  requiring  even  more  atomic 
data.  We  have  a  56  level,  steady-state  CRE  code  to  calculate  the  ionization 
dynamics  of  a  carbon  plasma.  This  code  gives  ion  fractional  abundances 
versus  temperature  with  the  total  ion  density  as  a  parameter.  See  section 
3.0  of  this  report  for  a  discussion  of  this  code,  its  applications  to  modeling 
our  experimental  results,  and  other  theoretical  work. 

ii)  Modeling 

A  population  inversion  can  only  exist  in  non-equilibrium  plasmas,  thus 
the  importance  of  assessing  the  plasma  state.  For  example,  population 
inversions  can  occur  in  a  supercooled  plasma,  partially  defined  by  the 
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property  that  >Te.  In  this  plasma  type  of  plasma  the  degree  of  ionization 
is  much  higher  than  that  implied  by  the  electron  kinetic  temperature. 
Inversions  can  occur  either  in  a  transient  ionization  build-up  phase  or  in  a 
recombination  phase.  Both  have  been  demonstrated  with  electron  beam 
pumping  but  in  the  latter  case  it  is  possible  that  a  steady  state  supercooled 
plasma  can  be  achieved.  Other  promising  non-equilibrium  situations  have 
been  demonstrated  in  the  distribution  of  bound  electronic  states.1-6 

We  have  performed  very  preliminary  initial  analysis  of  the  measured 
spectra.  Electronic  temperatures  have  been  obtained  from  Boltzmann  plots 
for  some  of  the  spectra.  These  are  based  on  a  rearrangement  of  the 

Boltzmann  equation  for  bound  electronic  state  densities  of  a  common 
ionization  state.  Substituting  Equation  1.14  into  Equation  1.11,  lumping  the 
constants,  and  taking  the  logarithim  of  both  sides  one  finds: 

ln(Ij^y/Ajjgj)  =  constant  -  Ej/kT  (1.19) 

A  plot  of  ln(IX/Ag)  versus  upper  state  energy  Ej  will  yield  a  straight  line  if 
those  excited  state  populations  follow  the  Boltzmann  equation.  An 
important  possible  source  of  error  in  this  analysis  is  the  assumption  of 

optically  thin  lines.  Other  errors  arise  due  to  uncertainties  in  A,  often  quite 
large.  In  order  to  obtain  meaningful  Boltzmann  plots  the  measured  spectra 
must  satisfy  the  following  conditions: 

i)  The  line  is  strictly  or  primarily  identified  with  one  atomic 

species  and  one  particular  transition  or  multiplet.  If  other 
components  or  transitions  are  coincident  or  close  (within 

instrumental  resolution)  ,  they  may  have  to  be  excluded  from 
the  calculation. 

ii)  The  signal  to  continuum  ratio  must  be  adequate  for  each  peak. 

iii)  The  transition  probability  must  be  available  for  each  line. 

One  line  can  often  be  associated  with  a  particular 
g-degenerate  multiplet,  which  cannot  be  separated  because  of 
instrumental  resolution. 

iv)  Electronic  temperature  estimates  from  line  ratios  require  at 
least  2  peaks  satisfying  requirements  (i)  -  (iii)  with 
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different  AEupper  (preferably  a  large  difference  to  reduce 
error).  Many  such  peaks  are  needed  to  get  an  adequate 
Boltzmann  plot. 

If  3  or  more  peaks  satisfying  these  requirements  are  available,  the 
negative  inverse  of  the  slope  of  the  linear  best  fit  line  gives  the  excitation 
temperature  for  that  species.  Electronic  equilibrium  for  that  species  exists 
if  the  correlation  coefficient  of  that  line  is  high  (>  0.9).  Because  of  the 
resolution  limits  of  the  current  0.275  meter  spectrograph  only  a  few  of  the 
spectra  obtained  satisfy  these  requirements.  We  expect  this  limitation  will 
be  removed  when  our  new  1.0  meter  spectrograph  becomes  available  in  the 
next  couple  of  months. 

If  the  plasma  components  all  have  similar  excitation  temperatures,  this 
suggests  that  this  temperature  may  also  be  the  electron  temperature.  For 
example  in  shot  M835  the  Boltzmann  analysis  yields  TexQu  =  2.4  ±  0.3  ev 
(Figure  1.11),  and  Tex0m=  2.7  ±  0.2  ev  (Figure  1.12).  Since  these  are 
identical  within  experimental  error,  our  initial  method  of  spectral  analysis 
assumed  that  the  electron  temperature  is  =  2.5  ev  and  that  the  plasma  is  in 
LTE.  The  use  of  the  Saha  equation  is  then  valid.  With  the  approximate  charge 
state  distribution  in  spectra  M835,  we  find  the  electron  density  = 
1018cnr3.  This  density  is  consistent  with  LTE  but  care  must  also  be  taken 
to  assure  that  the  argument  does  not  become  circular. 

The  assumption,  made  in  the  above  analysis,  that  Te  =  TJe  = 
TexCu  =  Texcni  is  strictly  unjustified.  Another  approximate  method  we 
have  used  to  estimate  the  plasma  state  uses  the  CRE  code  and  compares 
experimental  and  predicted  relative  ion  abundances.  The  method  involves 
finding  an  electron  temperature  required  to  give  a  charge  state  abundance 
consistent  with  experiment  for  a  particular  density.  This  procedure  is  also 
not  strictly  correct  as  the  presence  of  particular  lines  in  a  spectrum  does 
not  directly  translate  to  ion  abundances.  However  the  technique  can  yield 
some  useful  qualitative  information  about  the  plasma  state  and  plasma 
properties. 

The  spectra  for  shot  M835  and  M836  (taken  at  similar  times  in  the  pulse 
evolution)  indicate  that  CIII  and  CII  are  the  principle  constituents,  with 
enough  CIV  to  also  contribute  to  line  emission.  CIII  is  probably  more 
abundant  than  CII  as  the  CIII  lines  detected  have  high  upper  state  energies 
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(30-40  ev),  which  would  only  be  expected  in  a  hotter,  more  highly  ionized 
plasma.  Thus  we  are  looking  for  a  density  and  temperature  consistent  with 
abundances  decreasing  in  this  order:  CHI,  CII,  CIV.  We  expect  negligible  Cl 
and  CV  to  be  present.  Crossover  temperatures,  Tc,  can  be  defined  at  the 
points  where  various  ion  abundances  are  equal  (see  Figurel  .13).1  3  There  are 
two  crossover  temperatures  describing  the  region  where  a  particular  ion 
species  dominates.  For  example,  for  T  >  TCII:CII1  of  the  crossover  CII  to  CIII, 
the  ion  abundance  of  CIII  is  greater  than  that  of  CII.  Likewise  for  for  T  < 
T c 1 1 1 : c i v  tbe  crossover  of  CIII  to  CIV,  the  CIII  abundance  is  greater  than 
the  CIV  abundance.  Hence,  in  the  temperature  range  TCII.C1II  <  T  <  TC1IICIV, 
CIII  is  the  dominant  plasma  constituent.  Finally  for  T  <  TCII.CIV  for  the 
crossover  of  CII  to  CIV,  the  ion  abundance  of  CII  is  greater  than  that  for 
CIV.  These  three  temperatures  are  ordered  TC1I;CIII  <  TCII  CIV  <  TCIIICIV.  The 
temperature  range  required  in  order  to  be  consistent  with  experiment  is 
T c iirCiii  <  T  <  Tciiiciv*  Table  1.1  contains  these  crossover  temperatures, 
obtained  from  Figures  1.13  to  1.161-3,  for  four  different  densities. 

Table  1.1  shows  that  the  allowed  range  for  the  actual  electron 
temperature  (provided  by  the  first  two  columns)  is  larger  than  the 
excitation  temperatures  (2-3  ev)  for  these  four  densities.  Considering  the 
experimental  error,  the  case  for  Ne  *  1017  cm'3  is  close  to  agreement.  An 
upper  bound  on  the  temperature  is  provided  by  the  third  column,  the 
temperature  at  which  CIV  becomes  more  abundant  than  CIII.  These 
arguments  provide  some  evidence  that  equilibrium  assumptions  may  not 
apply  to  this  plasma,  and  the  actual  electron  temperature  may  be  somewhat 
higher  than  the  excitation  temperature. 

Table  1.1 
Tcrossover(ev) 

Ne(Cm-3)  onm I  CILCIV  CIILCIV 

1012  4  7  10 

1015  3  5  8 

1017  345 

1019  456 

A  better  way  to  utilize  the  CRE  code  is  to  first  obtain  an  independent 
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experimental  determination  of  electron  density,  perhaps  through  line 
broadening.  This  density  combined  with  a  relative  or  absolute  measurement 
of  line  intensities,  can  be  compared  with  the  CRE  prediction  for  those 
particular  lines,  to  determine  the  temperature.  Several  different  lines  or 
pairs  of  lines  can  be  used  to  cross  check  the  results.  In  addition  an 
independent  measurement  of  Te  would  be  desirable.  At  present  our 
instrumental  broadening  dominates  the  peak  widths  for  the  carbon  species. 
An  estimate  using  the  Hp  line  (hydrogen  being  the  principle  impurity),  may  be 
possible  if  it  can  be  separated  from  close  carbon  lines.  Stark  broadening  is 
expected  to  be  very  large  for  hydrogen  (and  hydrogen-like)  species.  These 
possibilities  are  expected  to  become  accesible  when  our  new  1.0  meter 
spectrograph  becomes  available. 
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1 .3)  Discussion  of  Spectroscopic  Results 


This  section  of  the  report  contains  complete  diode  analysis  and  spectra 
for  a  selection  of  crowbarred  and  noncrowbarred  shots.  The  diode  analysis 
consists  of  four  plots:  voltage,  current  with  current  model,  closure 
analysis,  and  dose,  just  as  in  Figure  1.7.  Recall  section  1.2.a.iv  and  the 
discussion  on  idealized  data.  The  OMA  gate  information  is  displayed  on  the 
voltage,  closure  and  dose  plots.  The  time  scale  for  the  OMA  is  referenced 
with  t  =  0.0  at  the  initial  rise  of  the  voltage.  Presented  on  the  closure 
plots  are  the  closure  velocity  and  the  predicted  diode  shorting  time.  The 
spectra  are  labeled  with  the  measured  wavelengths,  peak  identification, 
stage  identification,  and  excitation  temperatures  when  available.  All 
spectra  are  presented  on  the  same  relative  intensity  scale. 

a)  Crowbarred  Shots 

The  first  two  shots  M713  and  M718,  Figures  1.17  and  1.18,  are 
examples  of  stage  I  crowbarred  shots  where  no  light  is  detected.  For  both 
shots  the  dose  during  the  gate  is  below  200  J/gr.  This  is  well  below  the 
expected  lower  limit  for  anode  plasma  formation,  about  350  J/gr. 

In  shot  M717  the  light  detected  was  primarily  singly  ionized  carbon 
C 1 1 ,  with  some  molecular  emission  from  C2.  Note  the  low  level  of 
background  emission.  The  spectrum  is  shown  in  Figure  1.19a.  Tex,  the 
excitation  temperature  is  about  2.2  ev.  The  OMA  was  gated  from  52  to  631 
ns,  nominally  600  ns  wide.  The  diode  analysis,  in  Figure  1.19b,  shows  that 
the  dose  has  increased  to  >  370  J/gr  before  the  end  of  the  OMA  gate  so 
there  is  evidence  consistent  with  the  presence  of  an  anode  plasma.  The 
closure  analysis  implies  a  shorting  time  of  720  ns,  which  is  after  the  end 
of  the  gate  at  631  ns.  One  might  therefore  conclude  that  this  is  an  example 
of  a  stage  II  crowbar.  The  closure  velocity  analysis,  however  shows  that 
there  may  be  a  change  in  the  velocity  as  if  anode  plasma  closure  was 
coming  into  play.  The  current  model  with  bipolar  effects  does  not  matcK 
experiment  well  when  this  apparent  change  in  velocity  occurs.  There  are 
several  possibilities.  The  first  is  that  the  change  is  only  an  artifact  of 
experimental  error  although  there  still  is  an  anode  plasma.  There  may  not 
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Figure  1.17 
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in  fact  be  an  anode  plasma.  Another  possibility  is  that  a  simple  0-d  static, 
lumped,  uniform  plasma  model  probably  cannot  effectively  simulate  the 
diode-generator  and  plasma  dyanamics  with  anode  effects  in  3-d 
geometry.  A  more  complicated  current  model  would  be  required.  A  final 
possibility  is  that  this  is  a  transition  region  where  arc  and  space  charge 
effects  contribute  to  the  current.  If  this  does  indeed  represent  anode 
plasma  closure  then  the  diode  shorts  at  about  600  ns,  less  than  the  end  of 
the  OMA  gate.  Thus,  in  this  case,  and  because  of  our  uncertainty  in  the 
diode  dynamics  during  this  period,  the  light  may  also  be  contributed  from  a 
short  section  of  early  stage  III. 

Shot  M719  spectra  and  analysis  are  shown  in  Figures  1.20a  and  1.20b.  In 
this  shot  the  OMA  was  gated  from  461  ns  to  978  ns,  about  400  ns  later 
than  M717  with  a  nominally  500  ns  wide  gate.  Keep  in  mind  that  direct 
comparisons  between  shots  are  not  justifed  unless  the  closure  velocity 
and  effective  diode  area  are  close.  Due  to  the  irreproducability  of  the 
cathode  plasma  this  is  often  not  the  case.  The  shorting  time  in  this  shot  is 
about  the  same,  710  ns  which  places  it  in  about  the  middle  of  the  gate. 
Possible  contributions  to  the  detected  light  are  therefore  anode  plasma, 
cathode  plasma  and  arc  components. 

The  dose  from  space  charge  limited  current  incre  ases  to  >  320  J/gr. 
Recall  that  this  is  a  conservative  estimate  of  the  space  charge  dose  due  to 
use  of  a  simple  current  density  model  but  also  due  to  uncertainty  in 
assigning  a  cutoff  for  the  dose  integral.  This  cutoff  was  chosen  so  that  the 
current  model  agreed  with  experimental  current,  hence  space  charge 
assumptions  were  justified.  After  this  cutoff,  in  the  transition  and  arc 
region,  a  very  conservative  dose  estimate  is  used.  The  actual  dose  after 
this  cutoff  therefore  probably  increases  faster  than  shown.  Thus  there 
may  be  stage  II  anode  plasma  contributions  to  the  emitted  light.  Note  the 
presence  of  ionized  carbon  as  would  be  expected  from  an  anode  plasma 
with  incident  ionizing  electrons.  As  the  shorting  time  is  in  the  middle  of 
the  gate,  cathode  plasma  and  arc  components  would  also  be  visible  thus 
the  designation  of  early  stage  III.  The  molecular  components  are  all 
consistent  with  cool  hydrocarbon  cathode  contaminants.  The  increase  in 
maximum  detected  counts  (as  well  as  total  integrated  counts)  is  probably 
due  to  heating  by  the  high  current  (=10  kA)  low  voltage  arc  later  in  time. 

The  OMA  gate  for  shot  M720  is  from  976  ns  to  1494  ns,  approximately 
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500  ns  later  in  time  than  M719.  The  constituents  in  the  spectrum,  Figure 
1.21a,  are  identical  with  that  shot,  but  at  slightly  reduced  intensities. 
The  signal  to  noise  ratio  for  the  atomic  peaks  are  poor.  These  two  features 
are  consistent  with  a  plasma  that  is  cooling  and  recombining,  which 
correlates  with  the  falling  current.  The  diode  analysis  in  Figure  1.21b 
gives  a  shorting  time  of  902  ns,  hence  the  gate  is  after  the  arc  transition 
region,  and  hence  the  classification  of  middle  stage  III. 

The  last  three  shots  have  been  gated  sequentially  later  in  time.  Due  to 
the  width  of  the  gates,  they  make  up  a  long  time  integrated  temporal 
evolution  study,  and  give  a  good  idea  of  what  the  plasma  components  are. 
The  uncertainty  in  the  source  of  the  plasma  due  to  the  width  of  the  gates 
is  not  critical  since  we  desire  a  more  highly  ionized  plasma.  Another 
difficulty  in  deciding  what  stage  a  particular  shot  is  in,  aside  from  wide 
OMA  gates,  is  that  the  crowbars  end  just  as  the  dose  is  getting  up  into  the 
"interesting  regime"  where  anode  plasmas  are  expected.  In  noncrowbarred 
shots,  as  discussed  before,  the  doses  before  shorting  will  be  larger. 

b)  Noncrowbarred  Shots 

Due  to  the  voltage  opening  behavior  discussed  earlier,  the 
interpretation  of  noncrowbarred  shots  is  not  necessarily  easier  than 
crowbarred  ones.  Shot  M831  is  an  example  of  this.  The  spectrum,  Figure 
1.22a  shows  low  intensity  atomic  emission  from  CII  and  CIII.  We  have 
tentatively  associated  the  background  with  Cerenkov  emission  by  high 
energy  electrons  slowing  down  in  the  acrylic  window.  This  feature  is 
confusing  as  its  appearance  on  a  given  experimental  run  seems  to  be 
random,  however  it  does  not  appear  when  gates  are  during  the  low  energy 
portions  of  the  pulse.  The  OMA  gate  from  0  to  950  ns  goes  beyond  the 
apparent  shorting  time  of  873  ns  (see  Figure  1.22b).  The  dose  before 
shorting  has  reached  >  550  J/gr  well  into  the  range  for  anode  plasma 
initiation.  The  background  noted  on  the  spectrum  is  consistent  with  the 
formation  of  the  plasma  during  stage  II.  Due  to  the  width  of  the  gate  the 
line  emission  did  not  necessarily  occur  during  stage  II  and  there  may  be 
early  stage  III  elements  in  this  plasma.  The  uncertainty  in  the  plasma 
dynamics  because  of  the  voltage  opening  behavior  and  the  inability  to 
model  the  current  would  p"t  uncertainty  in  the  identification  of  the 
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plasma  source  even  had  the  gate  been  of  smaller  width.  The  next  shot  is  an 
example  of  this. 

The  OMA  gate  on  the  next  shot  is  the  shortest  presented  in  this  report, 
from  1075  to  1325,  only  250  ns  wide.  The  shorting  time  is  1400  ns  (Figure 
1.23b),  hence  the  gate  ends  before  predicted  shorting.  The  opening  behavior 
seems  to  be  turning  on  and  off  from  =  950  to  1250  ns,  thus  preventing  a 
clear  interpretation  of  the  spectra.  The  dose  reaches  >  940  J/gr.  The 
measured  spectrum,  Figure  1.23a,  shows  many  highly  ionized  carbon  peaks, 
including  a  probable  CIV  line.  Note  the  change  in  closure  velocity  from  5.0 
to  2.7  cm/ps  in  shots  M831  and  M839  respectively,  demonstrating  the 
irreproduceability  of  the  diode  load. 

The  next  six  shots  make  up  another  long  time  integrated  temporal 
evolution  study  of  the  3  ps  following  the  shot,  in  two  different 
wavelength  regions.  A  specific  discussion  of  the  analysis  is  not  necessary 
as  it  is  clear  that  these  are  all  stage  III  shots.  The  six  shots,  their  OMA 
gates,  wavelength  regions  and  figure  numbers  are  listed  in  Table  1.2 
below. 

Table  1 .2 

SbQl  QMA  gate  (p?)  Wavelength  (A)  Figure 
M835  1.0  -  1.95  3597-5066  1 ,24a, b 

M840  2.05  -  3.1  3597-5066  1 ,25a, b 

M842  3.0  -  4.075  3597-5066  1 ,26a, b 

M836  1 .05  -  1 .95  4602-6043  1 ,27a, b 

M837  1.95-  2.95  4602-6043  1.28a,b 

M838  3.05  -  4.05  4602-6043  1 .29a, b 

Each  succesive  shot  has  a  gate  of  nominally  1.0  ps  wide,  with  an 
additional  1.0  ps  of  delay.  The  shot  M831  gated  0.0  -  0.95  ps  is  an  example 
of  a  spectrum  taken  during  the  first  1.0  ps.  Below  are  presented  various 
comparisons  of  total  integrated  light  intensity  vs.  time.  We  will  ignore 
randomly  varying  window  attenuation,  and  differences  in  calibration  and 
light  sensitivity  for  different  grating  settings.  The  actual  difference  in 
detectable  or  emitted  light  in  those  wavelength  regions  may  be  obscured 
by  these  effects.  Differences  may  also  be  due  to  changes  in  diode 
parameters. 
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In  Table  1.3  are  the  total  number  of  counts  detected  in  the  given 
wavelength  regions  (normalizing  out  differences  in  OMA  gain),  and  the 
principal  detected  species.  Table  1.4  presents  the  same  data  normalized  to 
the  largest  integrated  intensity  overall  (on  the  left)  and  largest  integrated 
intensity  for  that  wavelength  region  (on  the  right).  The  normalization  on 
the  left  shows  that  the  most  intense  emission  is  immediately  following 
the  pulse  in  the  3597  -  5066  region  (shot  M835).  This  is  also  the  period 
with  the  most  highly  ionized  plasma.  The  normalization  on  the  right  shows 
that  the  time  decay  of  the  light  intensity  is  similar  in  each  range.  This  is 
significant  as  the  shots  were  obtained  through  differently  clouded 
windows.  Note  also  that  the  sequential  time  points  were  not  obtained 
completely  sequentially  in  the  shots^  which  should  minimize  the 
systematic  error. 

Table  1.5  presents  the  total  integrated  intensity  for  the  wavelength 
region  4602  -  5066  A,  the  region  common  to  both  studies.  The  numbers  in 
parenthesis  are  the  fraction  of  the  total  intensity  measured  that  was  in 
this  region.  Discounting  the  randomly  varying  window  attenuation,  and 
randomly  varying  electron  pulse/arc  characteristics,  the  only  difference 
in  total  light  detected  at  the  different  settings  would  be  due  to  the 
difference  in  grating  sensitivity.  The  numbers  in  the  third  column  are  the 
fraction  of  light  detected  by  the  275  setting  (4602-6043  A)  relative  to 
that  detected  by  the  225  setting  (3597-5066  A).  If  these  numbers  are 
fairly  constant,  the  above  two  effects  are  either  constant  or  not 
significant.  It  appears  that  the  two  shots  in  the  3  -  4  ps  period  are  not 
consistent.  This  is  probably  an  effect  due  to  differences  in  diode  behavior. 

Table  1 .3 


A  Reaion: 

3597-5066 

4602-6043 

principal 

Time 

0  - 1  ps 

Number  of  Counts 
32,967 

M83t 

species 

1  -  2ps 

2,085,640 

M835 

1,008,775 

M836 

CIII.CII.CIV 

2  -  3  ps 

868,208 

M840 

507,170 

M837 

cu.cm 

3  -  4  ps 

707,713 

M842 

189,379 

M838 

CI.CII 
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Table  1. 4 


0  - 1  ns  0.01/0.02 


1  -  2  jiS  1.00/1.00  0.48/1.0 

2  -  3  ps  0.41/0.41  0.24/0.50 

3  -  4  jiS  0.34/0.34  0.09/0.19 


Setting 

Total  Range: 

Time 

0  - 1  ps 

1  -  2ps 

2  -  3  jis 

3  -  4  ps 

On  one  shot  (M839)  we  gated  the  OMA  just  after  the  pulse  from  1.08  - 
1.32  jis.  This  gave  an  integrated  intensity  of  62,989  counts.  Compare  to 
shot  M836  which  had  a  gate  from  1  -  2  ps  (4  X  as  long),  but  gave  an 
integrated  intensity  of  1,008,775  counts  (16  X  as  many).  From  the  current 
traces  for  the  two  shots  we  note  that  the  current  is  still  either  rising  or 
steady  from  1.25  -  2  ps,  thus  plasma  production  and/or  heating  may  still 
be  increasing.  The  time  region  from  1-2  ps  is  quite  interesting  due  to  the 
intense  light  emission  and  high  ionization  states  detected  (up  to  CIV)  and 
requires  further  study. 

The  gates  used  in  this  study  are  very  wide,  and  a  lot  of  integration  over 
the  plasma  evolution  is  performed.  Tentatively  it  seems  as  though 
recombination  begins  as  the  arc  current  falls.  Rising  arc  current  seems  to 
drive  the  higher  ionization  states.  This  study  used  wide  gates  to  get  an 
initial  idea  of  the  species  present.  Use  of  shorter  gates  will  alleviate 
some  of  the  problems  associated  time  averaging  and  with  swamping  of 
details  by  intense  lines.  Excitation  temperatures  versus  time  are  shown  in 
Figure  1.30,  derived  from  the  analysis  in  section  1.2.b.ii.  More  weight 
should  be  given  to  the  data  points  higher  on  the  plot  as  they  are  obtained 
with  more  lines  in  the  Boltzmann  analysis.  The  excitation  temperature 


Table  1.5 
225 

(3597-5066) 


In  tens  it 


10,155  (0.31) 
876,943  (0.42) 
374,357  (0.43) 
411,247  (0.59) 


2Z5 

(4602-5066) 


636,659  (0.63) 
278,155  (0.55) 
111,151  (0.59) 
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Figure  1.30 


appears  to  be  constant  in  time  within  experimental  error. 

The  data  from  Figures  1.24b  to  1.29b  have  been  compiled  in  Table  1.6  to 
illustrate  the  irreproduceability  of  the  shots. 

Table  1.6 


Shot 

U  (cm/us) 

Reff  (cm) 

Tend  (ns) 

Tshort  (ns) 

TendATshoRT 

M831 

4.35 

2.91 

652 

874 

0.75 

M835 

3.96 

2.55 

803 

960 

0.84 

M836 

data  not  available 

M837 

3.32 

3.68 

635 

1145 

0.56 

M838 

3.99 

2.86 

456 

952 

0.48 

M839 

2.71 

3.60 

1089 

1402 

0.87 

M840 

5.3 

3.25 

610 

717 

0.85 

M842 

3.75 

2.77 

876 

1013 

0.87 

In  the  transition  region  between  stage  II  and  stage  III  we  believe  that 
there  is  both  arc-like  behavior  and  beam  behavior  (with  V  =  100  kV).  From 
the  above  analysis  this  is  the  region  with  the  most  intense  light  emission 
and  the  highest  ionization  states.  We  have  measured  up  to  CIV.  Higher 
ionization  states  can  be  detected  only  at  shorter  wavelengths  than  we 
have  explored.  We  intend  to  explore  this  transition  region  in  detail.  Shot 
M835  is  an  example  of  a  spectrum  taken  in  this  region. 

The  doses  during  the  initial  space  charge  regime  are  all  >  550  J/gr,  and 
keeping  in  mind  the  uncertainty  in  the  transition  region  and  arc  dose 
model,  the  doses  approach  2  kj/gr.  Although  the  mechanism  for  plasma 
formation  and  heating  in  the  arc  are  essentially  different,  reaching  these 
high  doses  during  the  space  charge  portion  of  the  pulse  could  be  very 
important.  As  discussed  above  and  in  the  conclusion  section,  there  are 
several  methods  we  can  use  to  do  this.  The  plasmas  created  in  stage  III 
are  also  interesting  themselves  as  tentative  analysis  (see  section  l.3.b.ii) 
indicates  they  may  be  nonequilibrium  plasmas.  The  high  currents  and  low 
energies  in  stage  III  result  in  substantial  coupling  giving  very  intense 
light  emission  from  a  highly  ionized  plasma.  This  is  probably  due  to  both 
heating  of  stage  II  plasma  inventory  and  the  creation  of  new  plasma 


inventory  through  anode  spot  mechanisms. 


1 .4)  Conclusions  on  Electron  Beam  Pumped  Carbon  Experiments 

In  stage  I  the  deposited  electron  dose  is  smaller  than  that  required  for 
anode  plasma  formation.  No  light  is  detected  in  this  phase. 

In  stage  II  the  deposited  electron  dose  is  large  enough  for  anode 
plasma  formation  before  diode  shorting.  Spectroscopic  data  for 
noncrowbarred  shots  (dose  =  500-800  J/g)  indicate  the  presence  of  CII 
and  CIII  with  a  high  level  of  background  radiation. 

In  stage  III  the  cathode  plasma  begins  shorting  the  anode  cathode  gap 
or  anode  and  cathode  plasmas  meet.  For  crowbarred  shots  (dose>300  J/g) 
moiecular  emission  from  C2  and  CH  and  atomic  emission  from  CII  and  H  is 
observed.  Noncrowbarred  shots  exhibit  intense  atomic  emission;  at  early 
times  CII  and  CIII  are  the  primary  components,  however  emission  from 
CIV  has  also  been  observed.  Excitation  temperatures  from  Boltzmann 
plots  are  2-3  eV.  At  later  times  the  emission  is  primarily  from  Cl  and  CII 
with  Tex  =  2-3eV.  The  electron  density  is  found  to  be  «  1018  cm'3  from 
the  Saha  equation. 

Future  experiments  will  emphasize  increasing  the  current  density  to 
obtain  higher  anode  doses.  The  principal  method  for  increasing  the  current 
density  will  be  the  use  of  an  external  magnetic  field.  Diode  aspects  which 
will  increase  the  current  density  are  decreasing  the  gap,  and  redesign  of 
the  cathode/anode  geometry. 


2.0)  Spectroscopy  of  Electron  Beam  Interactions  with  Helium 

Visible  emission  spectroscopy  was  performed  on  electron  beam 
produced  helium  plasmas  using  the  experimental  set-up  shown  in  Figure 
2.1.  The  Febetron  electron  beam  accelerator  was  used  for  these  intial 
studies  because  it  operates  with  a  higher  repetition  rate  and  does  not 
destroy  the  anode  foil,  permitting  many  more  shots. 

Typical  emission  spectra  for  helium  is  shown  in  Figure  2.2.  We  did 
not  observe  the  468.5  nm  (4-3)  transition,  which  had  been  observed  by 
Bunkin  et  al.2-1  to  undergo  stimulated  emission.  All  of  the  emitted 
spectral  lines  were  identified  as  neutral  helium.  The  electronic 
temperature  calculated  by  comparing  the  individual  helium  line  intensities 
was  found  to  be  =  0.8  ev.  This  is  much  lower  than  the  2  ev  electronic 
temperature  measured  for  argon  and  krypton  under  similar  experimental 
conditions.  Further  studies  are  in  progress  to  account  for  this  discrepancy 
and  also  to  investigate  the  production  of  Hell.  We  are  currently  running 
high  pressure  (100  torr  to  2  atms)  experiments  to  determine  the  effect  of 
pressure  on  temperature  and  degree  of  ionization. 

During  the  course  of  examining  the  emission  spectra  from  helium 
plasmas  produced  by  a  relativistic  electron  beam  we  found  that  if  very 
small  amounts  of  air  (as  small  as  50  millitorr  or  less)  were  present  as  an 
impurity  in  the  interaction  chamber  that  two  very  large  intense  lines  were 
observed  at  391.4  and  427.8  nm.  These  two  lines  correspond  to  the  (0-0) 
and  (0-1)  vibrational  transitions  in  the  first  negative  band  of  N2+,  B2iu  + 
-  x2Ig  +  .  When  the  N2  lines  were  present,  the  lines  originating  from 
helium  had  negligible  intensities!  We  pursued  this  important  effect  as  is 
described  in  the  following  paragraphs  and  published  in  Applied  Physics 
Letters,  (reprint  attached). 

We  very  carefully  added  measured  amounts  of  N2,  which  ranged  from 
0.1  to  10  Torr  partial  pressure  and  helium,  partial  pressures  from  0.1  to 
500  Torr.  In  time-resolved  studies  performed  with  the  Optical 
Multichannel  Analyzer  (OMA),  gates  of  50  ns  delayed  by  20  ns  for  each 
pulse,  the  N2+  light  emission  was  found  to  follow  the  beam  voltage  after 
the  initial  rise  time  (see  Figure  2.3).  During  the  time  resolved  studies,  no 
other  lines  were  observed  except  for  the  391.4  nm  and  427.8  nm  nitrogen 
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bands.  Figure  2.4  shows  the  relative  intensity  of  the  391.4  nm  band  as  a 
function  of  helium  to  nitrogen  ratio  for  nitrogen  partial  pressures  of  0.1,  1 
and  10  Torr.  Figure  2.5  shows  similar  results  for  the  427.8  nm  band. 
Figure  2.6  shows  emission  spectra  for  a  variety  of  helium  pressures. 

As  seen  in  Figs.  2.4  and  2.5,  the  largest  light  intensities  at  391.4  and 
427.8  nm  for  a  given  pressure  of  nitrogen  occur  for  the  largest  ratio  of 
He/N2,  i  e.  the  smallest  amount  of  N2.  Also,  the  more  helium  present,  the 
more  light  is  produced  up  to  total  pressures  of  about  300  Torr  where  the 
intensity  levels  off.  The  linear  portion  of  the  curves  correspond  to  a  slope 
of  about  1  for  the  range  of  total  pressure  10  -  300  Torr  for  the  cases 

where  the  N2  partial  pressure  was  1  and  0.1  Torr.  The  slope  was  0.5  when 
the  partial  pressure  of  N2  was  10  Torr.  For  cases  where  the  helium 
partial  pressure  is  smaller  or  about  the  same  as  the  nitrogen  partial 
pressure,  there  is  very  little  increase  in  the  N2+  lines  compared  to  the 
case  of  nitrogen  alone. 

As  seen  in  Fig.  2.2,  the  largest  intensity  lines  observed  when  no  added 
nitrogen  was  present  were  388.8,  447.1,  501.5,  and  587.5  nm.  In  the  case 
of  the  388.8  nm  line,  the  energy  of  this  transition  is  very  close  to  that  of 
the  N2+  391.4  nm  band.  Therefore,  a  possible  mechanism  for  the 
excitation  of  this  band  is 

N2+  (X2I)  +  He  (33P)  -  >  N2+  (B2!)  +  He  (23P). 

Since  the  slope  of  the  log  of  the  intensity  as  a  function  of  the  log  of  He/N2 
is  close  to  1  for  small  N2  partial  pressure  cases,  there  appears  to  be  a  one 
to  one  correspondence  between  the  number  of  photons  emitted  for  the 
given  bands  and  the  number  of  helium  atoms  present,  above  a  particular 
helium  to  nitrogen  threshold.  This  one  to  one  correspondence  is  cc  isistent 
with  the  above  reaction  and  a  simple  kinetics  model.  The  details  of  the 
kinetics  model  are  given  in  the  attached  reprint. 

For  formulating  a  model  of  four  simple  reactions,  we  can  solve  for 
the  time  rate  of  change  of  the  He(3P)  state,  which  is  then  used  to  find  the 
time  rate  of  change  of  N2+(B2i).  From  these  solutions  we  can  show  that 
the  emission  rate  for  N2+  (B2!  -  X2I)  is  linearly  proportional  to  the 

electron  beam  density  (which  we  observe),  the  ground  state  population  of 


N  2+  and  most  importantly  linearly  proportional  to  the  He  ground  state 
concentration.  This  is  an  important  result  because  other  investigators  who 
have  examined  the  effects  of  helium  on  the  N2+  laser  have  felt  that  the 
N2+  emission  is  the  result  of  interactions  with  He2  or  He2+  or  more 
complicated  processes  which  turn  out  to  be  precluded  by  our  results. 

We  recently  published  these  results  in  Applied  Physics  Letters  and 
presented  a  paper  at  the  38th  Annual  Conference  on  Gaseous  Electronics 
where  they  were  well  received. 

As  mentioned  previously,  now  that  we  have  fully  investigated  the 
effects  of  air  impurities  on  our  helium  emission  spectra,  we  are  setting 
up  our  experiment  to  eliminate  this  effect  so  that  we  can  examine  the 
interaction  of  the  electron  beam  with  helium  alone.  We  plan  to  investigate 
the  possibility  of  stimulated  emission  in  helium,  particularly  for  high 
pressures  (>  1  atm). 


Detector 


Figure  2.2 


3.0)  Theoretical  Research 


This  is  a  summary  of  the  work  being  done,  as  of  November  1986,  to  explain  the 
experimental  observations  of  a  carbon  plasma  heated  by  an  incident  electron  beam.  An 
overveiw  of  the  theory  is  presented  as  well  as  a  description  of  proposed  work. 

The  main  goal  of  the  investigation  is  to  formulate  a  theoretical  model  which  will 
describe  the  motion  of  the  carbon  plasma  and  predict  the  spectrum  of  emitted  light.  The 
predicted  spectra  may  then  be  directly  compared  to  spectra  which  are  experimentally 
measured.  The  model  being  currently  developed  is  based  on  work  done  by  J.W. 
Thornhill.3-1  Three  main  areas  must  be  addressed:  hydrodynamic  description  of  the 
plasma,  an  ionization  dynamics  model,  and  the  energy  deposition  of  the  incident  particle 
beam.  The  work  done  in  each  of  these  areas  is  presented  and  explained.  Due  to  the 
complexity  of  the  problem  numerical  methods  are  used  to  obtain  results.  The  code 
which  is  being  developed  is  not  described  here. 

The  main  hydrodynamic  assumption  is  that  the  plasma  can  be  described  using 
single  fluid  and  temperature  equations.  This  places  a  constraint  on  the  time  scale  of  the 
problem.  Namely  that  the  electron-ion  thermalization  time  should  always  be  smaller 
than  the  hydrodynamic  time.  This  constraint  is  not  serious  but  should  be  kept  in  mind. 
The  model  is  developed  in  one  dimensional  lagrangian  coordinates.  This  simplifies  the 
conservation  equations  and  allows  for  simulation  of  situations  where  extensive 
compression  or  expansion  occurs. 

The  equation  of  motion  (Eq.  3.1)  is  obtained  by  performing  a  force  balance  on  an 
incremental  volume  of  plasma.  The  plasma  is  assumed  to  be  inviscid.  The  Von 
Neumann  artifical  viscosity  is  introduced  to  allow  finite  difference  techiques  to  be  used  in 
situations  where  discontinuities  in  pressure  and  density  occur.3-2 


dll  sj  d  1 

_  =  x  -  P  +  q  + - AMOM 

at  am  L  J  p5v 


(3.1) 


u  -  fluid  velocity  q  -  Von  Neumann  artificial  viscosity 

P  -  plasma  pressure  5v  -  volume  of  incremental  element 

p  -  plasma  density  m  -  Lagrangian  mass 

AMOM  -  momentum  imparted  by  beam  electrons  to  plasma 
5  -  1  for  cartesian  geometry,  2  for  cylindrical,  3  for  spherical 
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The  energy  conservation  equation  (Eq.  3.2)  is  obtained  by  performing  an  energy 
balance  on  the  incremental  volume  element.  It  is  assumed  that  Fourier's  Law  accurately 
describes  the  heat  transfer.  No  flux  limiter  is  needed  since  the  electron  mean  free  path  is 
assumed  smaller  than  the  temperature  gradient.  This  should  be  true  in  the  plasmas  we 
are  considering.  The  final  form  after  some  manipulation  is  given  below. 
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T  -  plasma  temperature  V  -  specific  volume 

S  -  energy  source  to  plasma  K  -  thermal  conductivity 

Cv  -  specifc  heat  at  constant  volume  ( )v  -  implies  at  constant  volume 
R  -  radiation  emitted/absorbed  by  plasma 


The  ionization  dynamics  model  is  used  to  calculate  the  equation  of  state 
information  that  is  needed  in  the  hydrodynamic  simulation.  The  parameters  of  interest 
are  the  internal  energy,  the  effective  charge,  and  the  radiated  power  from  free-free, 
free-bound,  and  bound-bound  transitions.  This  information  is  obtained  in  tabular  form  as 
a  function  of  temperature  and  ion  density  and  supplied  as  input  to  the  hydrodynamics 
model.  In  addition,  the  ionization  dynamics  model  is  solved  concurrently  with  the 
hydrodynamics  model  to  provide  current  state  populations  and  the  plasma  spectrum. 
This  arrangement  allows  the  calculation  of  both  spatially  and  temporally  dependent 
spectra. 

The  Collisonal-Radiative  Equilibrium^  (CRE)  model  is  used  to  describe  the 
ionization  dynamics  of  the  system.  This  model  consists  of  a  set  of  rate  equations,  one  for 
each  quantum  state  included,  which  describe  the  effect  of  atomic  processes  on  the 
population  and  depopulation  of  a  given  state.  The  model  used  is  an  equilibrium  model 
since  no  time  dependence  is  included  in  the  rate  equations. 

Our  version  of  the  CRE  model  has  a  few  additional  restrictions  placed  on  it.  The 
effects  of  photon  stimulated  emission  and  absorption  are  not  included.  This  is  done  to 
avoid  the  complication  of  including  radiation  transport  and  having  to  solve  a  nonlinear 
set  of  equations.  The  nonlinearity  is  a  result  of  the  direct  coupling  of  the  radiation 
transport  and  rate  equations.  The  neglect  of  these  effects  restricts  the  applicability  of  the 
model  to  optically  thin  plasmas. 

Lowering  of  the  ionization  potential  is  not  explicitly  included.  This  phenomenon 
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sets  a  bound  on  the  number  of  discrete  quantum  states  available.  It  implicitly  affects  the 
model  by  allowing  use  of  a  finite  number  of  states  be  adequate  to  describe  the  plasma. 

The  atomic  processes  included  in  our  model  are:  collisional  ionization,  three  body 
recombination,  collisional  excitation,  collisional  de-excitation,  spontaneuos  emission, 
radiative  recombination,  and  dielectric  recombination.  The  rate  coefficients  are 
calculated  using  various  methods  and  the  ones  we  are  using  were  obtained  by  J.W. 
Thornhill  from  D.  Duston  at  the  Naval  Research  Laboratory  in  Washington,  D.C.  Recent 
work  by  D.  Duston  summarizes  how  the  rate  coefficients  were  calculated.3-4 

The  model  currently  includes  56  states.  It  should  be  noted  that  many  of  these 
states  actually  represent  several  degenerate  levels  of  carbon  which  have  been 
averaged  over  the  angular  momentum  and  magnetic  quantum  numbers.  The  population 
of  the  degenerate  levels  can  be  recaptured  by  assuming  that  they  are  in  equilibrium  and 
using  a  partition  function  to  extract  their  populations  from  the  CRE  calculated  state 
population. 

The  final  component  of  our  theory  is  the  energy  deposition  model.  The  approach  is 
simplistic  at  the  present.  The  stopping  power  of  an  incident  relativistic  electron  is 
integrated  over  the  plasma.  This  is  then  multiplied  by  the  current  density  of  the  beam  at 
each  point  in  the  plasma  to  obtain  the  desired  result.  The  stopping  power  is  based  on 
the  work  by  D.  Mosher  for  relativistic  electrons  interacting  with  a  plasma  which  has  both 
free  and  bound  electrons.3-5  A  collision  term  was  found  using  Fokker-Planck  formalism. 
Assuming  the  incident  beam  can  be  treated  as  a  delta  function  in  momentum  space,  the 
appropriate  moment  of  the  collison  term  yields  the  stopping  power  (Eq.  3.3).  Note  that 
the  energy  loss  due  to  electron-ion  collisions  and  bremsstrahlung  has  been  ignored. 


dE 

dx 


=  -47rNr2  ae  E0 


e2-E2 

o 


(3.3) 


ae  =  ZAln  (P2y(y-1)i/2E0  /  lz]  +  Zln[  XD\Z  /  2ytE0  ] 
P  =  v/c  Y=E/E0 


E  -  incident  electron  total  energy 
E0  -  electron  rest  mass  energy 
lz  ■  ionization  energy 
ZA  ion  atomic  number 


N  -  plasma  density 
r  -  classical  electron  radius 
XD  -  debey  length 
Z  -  effective  charge  of  plasma 


sop  °g  power  was  found  to  yield  ranges  which  are  in  good  agreement  with 
.'  standard  formulas.  Thus  it  is  assumed  to  be  accurate  enough  for  the 
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level  of  developement  thus  far.  It  should  be  noted  that  the  attenuation  of  the  beam 
electrons  by  the  plasma  is  not  accounted  for  when  calculating  the  energy  deoosition. 

Much  work  is  planned  and  under  way  to  revise  the  theory  and  hopefully  extend  its 
validity.  The  hydrodynamic  portion  of  the  code  is  being  reveiwed  and  rewritten  so  that  it 
is  more  applicable  to  describing  the  experimental  configuration.  Also  subroutines  which 
calculate  and  output  the  line  emission  spectrum  are  being  developed. 

In  addition  the  energy  deposition  model  will  be  reveiwed.  The  attenuation  of  the 
incident  electrons  as  they  tranverse  the  plasma  must  be  accounted  for  to  have  a 
satisifying  theory.  It  is  felt  that  an  exponential  attenuation  might  be  adequate  and 
consistent  with  the  level  of  the  model’s  sophistication. 

Another  problem  arises  when  one  attempts  to  compare  with  experiment.  It  is  not 
clear  how  to  calculate  the  incident  beam's  current  density  after  a  period  of  time.  This 
information  is  vital  to  the  energy  deposition  and  the  question  needs  to  be  addressed.  If  it 
becomes  possible  to  measure  the  current  density  instead  of  the  current,  then  the  current 
density  could  be  input  directly  to  the  code  and  the  problem  circumvented.  Also  at  this 
time  it  is  felt  that  most  of  the  energy  deposition  occurs  when  the  plasma  arcs  the 
anode-cathode  gap.  How  to  accurately  describe  the  energy  deposition  during  this 
phase  of  the  experiment  needs  to  be  considered  for  accurate  simulation  of  experimenta1 
conditions. 

The  ionization  dynamics  model  needs  to  be  carefully  reveiwed  The  CRE  model  is 
valid  for  dynamic  plasmas  with  free  electron  densities  above  1O10 ;cm3  This  threshold  is 
obtained  by  requiring  the  relaxation  times  for  atomic  processes  be  shorter  than  the 
hydrodynamic  time  scale  The  longest  relaxation  times  are  those  to  ground  states  and 
these  set  the  density  limit  It  desired  to  describe  plasmas  with  tree  electron  densities 
below  this  threshold  The  Quasi-Stationary  Collisionai  Radiative  Model  which  retains 
time  dependence  in  the  ground  state  rate  equations  might  be  used  to  extend  to  iower 
density  plasmas  ^  ®  This  is  one  model  which  will  be  considered  Note  that  if  the  plasma 
is  static  the  CRE  model  would  be  valid  at  iower  densities 

The  adequacy  of  the  number  of  states  used  -n  the  model  must  aiso  be  adressed 
One  method  to  include  additional  states  is  the  brute  tome  approach  o*  s>mp;y  addmq 
them.  This  approach  increases  the  number  of  rate  equat  ons  wnich  pe  snveo  a"d 
may  become  prohibitively  expensive  *rom  a  computahona  vcwpomt  Ancthe'  poss.b ■*> 

a  •* 

is  to  use  a  Partial  Local  Thermodynamic  Eq^  ;ipr  um  Mode  J  s  approach  uses  "  e 
reduced  ionization  potential  to  find  a  thresho  d  energy  above  w’  k.'1  the  quantum  states 
are  in  local  thermodynamic  equil  bnum  l  T F  One  uses  {  Tf  —ode1  ab've 
threshold  to  determine  the  state  population  .  Another  model  suid  as  tk  e  C.Rt  s  usee 


below  the  threshold  to  obtain  populations  of  the  low  lying  states. 

These  additions  and  alterations  to  the  current  theory  will  be  reveiwed.  The 
appropriate  changes  will  be  implemented  where  necessary. 
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Relativistic  electron  beam  interactions  with  very  small  ratios  of  nitrogen  to  helium  ( 10' 

10  4)  have  been  found  to  produce  extremely  large  N,*  {B'lf  -X  'If  )  intensities  at  391.4 
and  427.8  nm.  compared  to  line  intensities  originating  from  helium  These  results  occurred  in 
the  total  pressure  regime  of  0.1-500  Torr  The  pressure  scaling  results  presented  here  are 
inconsistent  with  previously  proposed  kinetic  mechanisms  for  tht  N,'  laser  pumped  by 
helium  With  a  simple  model  of  the  chemical  kinetics,  we  show  that  this  effect  is  due  to  the 
collisional  transfer  of  energy  between  excited  states  of  helium  atoms  and  the  ground  state  of 

N/- 


It  is  well  known  that  population  inversions  of  the 
B:1f  -X:1f  electronic  transition  of  NY  can  be  achieved 
in  the  presence  of  helium  1  Two  possible  mechanisms  have 
been  used  to  describe  the  pumping  process  Collins'  believes 
the  upper  laser  level  is  pumped  by  charge  transfer  from  He,' 
to  N,.  Kokawaera/.:  believe  that  electron  impact  excitation 
is  the  dominant  pump  mechanism  Most  of  the  research  per¬ 
formed  on  the  NY  laser  has  involved  large  overall  He-N, 
discharge  pressures  (1-36  atm),  :  where  the  presence  of 
He,'  may  be  feasible  Investigations  of  lower  (  <  50  Torr) 
overall  pressures,  however,  have  observed  population  inver¬ 
sions  in  atomic  nitrogen'4  rather  than  in  molecular  nitro¬ 
gen  Population  mechanisms  proposed  for  atomic  nitrogen 
lasers  range  from  excitation  of  nitrogen  atoms  by  collisions 
between  metastable  excited  atoms  of  helium  (2'S)  and  ni¬ 
trogen  molecules'  to  the  recombination  of  electron-ion 
pairs  4  In  this  letter,  we  present  intensity  measurements  for  a 
wide  range  of  helium  to  nitrogen  ratios  for  relativistic  elec¬ 
tron  beam  pumped  gases  in  the  intermediate  pressure  range 
ofO  1-500  Torr  and  show  that  these  results  are  inconsistent 
with  the  above  mentioned  mechanisms  We  also  present  a 
simple  chemical  kinetic  model  which  is  consistent  with  our 
experimental  observations. 

The  experimental  apparatus  has  been  described  pre¬ 
viously  '  The  experiment  consisted  of  a  relativistic  electron 
beam  at  a  peak  voltage  of  300  keV,  a  peak  current  of  1  kA. 
and  a  pulse  width  of  300  ns  The  diode  consisted  of  a  2  5-cm- 
diam  carbon  brush  cathode  and  a  0  025-mm  titanium  anode 
foil  The  beam  was  injected  into  an  aluminum  vacuum  vessel 
with  lucite  windows  on  the  sides  and  end  The  emitted  light 
at  several  helium  partial  pressures  each  with  nitrogen  partial 
pressures  of  0  1,  1,  and  10  Torr  was  sampled  axially  The 
light  was  analyzed  by  a  0  275-m  spectrograph  coupled  to  an 
optical  multichannel  analyzer  (OMA)  The  OMA  was  gat¬ 
ed  with  pulses  of  50  and  500  ns 

For  all  experiments  where  helium  was  present  with 
small  to  large  amounts  of  nitrogen,  two  very  intense  lines 
were  observed  at  391  4  and  427  8  nm  These  lines  correspond 
to  the  (0-0)  and  (0-1 )  vibrational  transitions  in  the  first 


negative  band  of  N2*  ,  B2lf-X2lf  Lines  originating 
from  atomic  or  molecular  helium  had  negligible  intensities 
for  these  cases.  When  no  nitrogen  was  present,  the  major 
lines  were  attributed  to  neutral  atomic  helium  at  388.8, 
447  1,  501.5,  and  587.5  nm  In  time-resolved  studies  per¬ 
formed  with  OMA  gates  of  50  ns  delayed  by  20  ns  for  each 
pulse,  the  N24  light  emission  was  found  to  follow  the  beam 
voltage  after  the  initial  nse  time  (see  Fig.  1 ).  During  the 
time-resolved  studies,  no  other  lines  were  observed  except 
for  the  391.4-  and  427.8-nm  nitrogen  lines  Figure  2  shows 
the  relative  intensity  of  the  391.4-nm  band  as  a  function  of 
helium  to  nitrogen  ratio  for  nitrogen  partial  pressures  of  0  1 . 
1,  and  10  Torr  Figure  3  shows  similar  results  for  the  427  8 
nm  band. 

As  seen  in  Figs  2  and  3,  the  largest  light  intensities  at 
391  4  and  427.8  nm  for  a  given  pressure  of  nitrogen  occur  for 
the  largest  ratio  of  He/N  ,  That  is,  the  more  helium  present, 
the  more  light  is  produced  up  to  total  pressures  of  about  300 
Torr  where  the  intensity  levels  off  The  linear  portion  of  the 
curves  corresponds  to  a  slope  of  about  1  for  the  range  of  total 
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HO  I  Beam  voltage  >•  \  and  W|  4  nm  (O)  and  *2~  H  nm  (•)  emission 
band  inienvmev  av  a  function  of  lime  The  partial  prevvurev  of  N  and  helium 
were  OS  and  100  Torr  Optical  multichannel  analwer  gates  were  X0n\  and 
nmev  pinned  are  for  ihe  cenier  of  the  gale 
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FIG  2  Inlensiiy  of  the  1  4-nm  emission  tend  as  a  funciior  of  He  ""N 
ratio  for  N,  pressures  of  0  1  Ton  (  «  ).  I  Ton  (O).and  lOTon  (Cl  The 
slopes  of  the  linear  ponions  of  (he  curses  are  as  indicated 

pressure  10-300  Torr  for  the  cases  where  the  N\  partial  pres¬ 
sure  was  1  and  0  1  Ton  The  slope  was  0.5  when  the  partial 
pressure  of  N,  was  10  Torr  For  cases  where  the  helium 
partial  pressure  is  smaller  or  about  the  same  as  the  nitrogen 
partial  pressure,  there  is  very  little  increase  in  the  intensity  of 
the  N,*  lines  compared  to  the  case  of  nitrogen  alone 

As  mentioned  previously,  the  largest  intensity  lines  ob¬ 
served  when  no  added  nitrogen  was  present  were  388  8. 
447  1.  501  5,  and  587  5  nm  In  the  case  of  the  388  8-nm  line, 
the  energy  of  this  transition  is  very  close  to  that  of  the  N  ‘ 
391  4  nm  band  Therefore,  a  possible  mechanism  for  the  ex 
citation  of  this  band  is 

NV  (T  ;1  l  *  He<  37*)  AV  iB:l)  -  Het  2  '/*).  tii 

w  here  k  ,  is  the  reaction  rate  constant  Since  the  slope  of  the 
log  of  the  intensity  as  a  function  of  the  log  of  He  N  is  close 
to  I  for  the  small  N  partial  pressure  cases,  there  appears  to 
be  a  one  to  one  correspondence  between  the  number  of  pho 
tons  emitted  for  the  given  bands  and  the  number  of  helium 
atoms  present,  above  a  particular  helium  to  nurogen  thresh 
old  This  one  to  one  correspondence  is  consistent  with  reas 


lion  ( 1 )  as  illustrated  by  a  simple  kinetic  model  Assume 
that  the  electron  beam  excites  the  helium  from  the  ground 
state  to  the  3  *5  state  as  follows 

He(  l‘S)  -f  •He(3'/>)  +  e,  (2) 

where  k:  is  the  rate  constant  for  reaction  (2 )  and  is  a  func¬ 
tion  of  the  electron  beam  energy  At  low  pressures,  the  main 
deactivation  mechanism  of  HeO'/*)  and  N}'  (B:l)  is  ra¬ 
diative  decay 

He(37»)-He(2'S)  4  AvliiilMll.  (3) 

N,*  (£  J2  (T :I)  +  Av„,  4 (4) 

Note  that  He(2'5)  is  a  metastable  state.  We  also  make  the 
following  assumptions  all  back  reactions  are  negligible, 
collisiona!  decay  of  He(3'/>)  is  unimportant  because  the 
exchange  reaction  He( 3'/*)  4  He(2'S)-He(2\S) 
4  He(  3'/*)  does  not  alter  the  He(3'P)  concentration,  the 
collisional  quenching  of  N/  [B'l)  is  unimportant  for  low 
pressures  (although  this  assumption  is  discussed  later),  and 
that  the  electron  beam  density  can  be  approximated  by  an 
exponential  of  the  form  [e„  ]  =  |r„  )"(<■  "  '  -  e  "') 
where  [c,.  ]  "is  some  constant  characteristic  of  the  beam  den 
sity  and  u1,  and  u  ■  are  constants  which  describe  ihe  beam 
current  tempioral  profile  W>  now  write  a  differential  equa 
tion  for  He(  3 'P)  as 

d  |Hri  T  ) 
d\ 

k  .(  Het  I  '5  )  ]  |r,  |"(r  “  .  - 

A  ,  |  Hct  T  7’ i  I  A  ,  |  Hr  i  V/’i  |  |  N  ,  \  1  ,  )  ,  «■  , 

Now  assume  that  ibr  ground  Male  populationv  Hr  i  I'S)  and 
N  .‘  (  X  2  )  stay  approximairls  constant  Ihrn  a  linear  ho 
mogenrous.  ordinary  diffrrrniial  equation  exists  foi 
|  He i  T/’l)  and  can  hr  solved  using  siandard  tnhniques 
The  solution  is 
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FIG  '  Intrmit*  of  the  42?  I  nm  band  rmiunm  as  a  fumitory  <4  Mr  *s 
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vx  here  A  A  ,  |  N  ■  I  i  >  )  «  A ,  and  4  A  {Me  I 
•  If,  |  We  tan  alv-  write  the  ratr  equation 
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A  rale  equation  for  the  time  rate  of  change  of  the  photon 
emission  at  391.4  nm  can  be  written  as 

••nl  =&,[N/  [B: 1)].  (9) 

dt 

Substituting  the  results  of  Eq  (8)  gives 

d  |  =  *,*«*,{  NV  {X:2 )  ]  ( He(  I '$) )  [e* )" 

dt  *  J 

*  (  («■;-«•,  He  k'-e  k-') 
l  (k  -  u\)(A  -  u’,)(A4  -  k) 


(A  -  w,)(kt  -  w, ) 
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Thus,  the  emission  rate  for  N/  ( B:l-X  :1)  is  linearly  pro¬ 
portional  to  the  He  ground-state  concentration,  the  electron 
beam  density,  and  the  ground-state  population  of  N/  If 
collisional  quenching  of  N,‘  (B  2 )  is  included  as 

k4 

N.-  (B'l)  -  M  — N/  (T'2)  -  M.  (ID 

then  A,in  Eq  I  10 )  would  be  replaced  by  k  '  =  k4  4  k,  \M  j, 
where  M  is  anv  collisional  panner  such  as  He  or  N,  Equa¬ 
tion  (10 1  would  then  have  a  mo  e  complicated  dependence 
upon  helium  concentration  than  the  simple  linear  one  shown 
here  Obviously,  to  besi  investigate  the  effects  of  collisional 
de-excitation  and  to  remove  ihe  assumption  that  N/  ( T  "2 ) 
and  He i  I  ,S  )  remain  constant  would  require  accurate  rate 
constants  and  a  computer  vdution  of  the  differential  equa 
lions 

Another  important  result  of  the  data  shown  in  figs  2 
and  t  is  that  as  the  nitrogen  partial  pressure  increases,  the 
overall  amount  of  line  enhancement  decreases,  as  well  as  the 
slope  of  the  linear  part  of  the  intensity  versus  pressure  ratio 
t  urvrs  for  nitrogen  partial  pressures  of  0  I  Torr.  the  emif 
led  iniensiis  increases  bv  almost  a  factor  of  100  whereas  for 
nitrogen  partial  pressures  of  10  Torr.  the  overall  increase  is 
less  than  a  factor  of  10  compared  to  the  no  helium  case  It 
appears  that  a  vers  important  collisional  depopulating 
mechanism  of  N  Iff  li  is  collisions  with  other  nitrogen 
rrtolec  ules 

Never  at  anv  lime  wav  emission  due  to  He  -  or  prrstulated 
Me  iserKel  h.  observed  l  his,  lit  addit  <>n  to  the  fact  that 
the  highest  .  cnail  pressme  examined  was  Vg)  Tori,  leads  to 


the  conclusion  that  charge  transfer  between  He/  was  prob¬ 
ably  negligible,  also  since  for  example  the  388.8-nm  helium 
line  which  populates  the  metastable  state  of  helium  was  not 
observed  to  increase  or  decrease  with  the  39 1  4-nm  band,  in 
contrast  to  the  results  of  Atkinson  and  Sanders. '  Instead,  the 
388.8-nm  helium  line  was  observed  only  in  the  absence  of 
nitrogen,  thus  it  is  doubtful  that  direct  transfer  between  the 
metastable  state  of  helium  (2:5')  and  the  N/  ( B  }1 )  was  sig¬ 
nificant.  Kokowa  et  a!.  -  claim  that  electron  impact  ioniza¬ 
tion  is  the  major  pumping  process  and  that  the  presence  of 
helium  only  helps  to  raise  the  electron  temperature  rather 
than  supplying  reactants  for  the  pumping  process.  The  re¬ 
sults  presented  here  show  a  definite  helium  dependency,  par¬ 
ticularly  since  we  have  shown  previously’  that  the  electron 
temperature  is  pressure  insensitive  for  our  e-beam  produced 
discharges 

Relativistic  electron  beam  produced  plasmas  in  helium- 
nitrogen  mixtures  produce  very  large  intensity  lines  at  391.4 
and  427.8  nm.  These  Lnes  originate  from  the  vibrational 
transition  of  (0-0)  and  (0-1 )  in  the  first  negative  system  of 
N/  .  The  line  intensities  increase  as  the  ratio  of  helium  to 
nitrogen  increases  The  greatest  increase  is  seen  for  the  low¬ 
est  nitrogen  concentrations  The  increased  intensity  is  linear 
with  helium  pressure  between  the  total  pressure  of  10-300 
Torr  for  small  partial  pressures  (0  l-1.0Torr)  of  nitrogen 
This  increase  is  believed  to  be  due  to  collisional  transfer  of 
energy  between  helium  excited  states  and  the  N/  ground 
stale  as  show  n  by  the  agreement  between  our  simple  kinetic 
model  and  experimental  results 
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